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ABSTRACT 



We study the formation of molecular precursors to dust in the ejecta of Population III 
supernovae using a chemical kinetic approach to follow the evolution of small dust 
cluster abundances from day 100 to day 1000 after explosion. Our work focuses on 
zero-metallicity 20 Mq and 170 M progenitors, and we consider fully-macroscopically 
mixed and unmixed ejecta. The dust precursors comprise molecular chains, rings and 
small clusters of chemical composition relevant to the initial elemental composition of 
the ejecta under study. The nucleation stage for small silica, metal oxides and sulphides, 
pure metal, and carbon clusters is described with a new chemical reaction network highly 
relevant to the kinetic description of dust formation in hot circumstellar environments. 
We consider the effect of the pressure dependence of critical nucleation rates, and test 
the impact of microscopically-mixed He + on carbon dust formation. Two cases of metal 
depletion on silica clusters (full and no depletion) are considered to derive upper limits to 
the amounts of dust produced in SN ejecta at 1000 days, while the chemical composition 
of clusters gives a prescription for the type of dust formed in Pop. Ill supernovae. 

We show that the cluster mass produced in the fully-mixed ejecta of a 17OM pro- 
genitor is ~ 25 M & whereas its 20 Mq counterpart forms ~ 0.16 Mq of clusters. The 
unmixed ejecta of a 170 Mq progenitor supernova synthesizes ~ 5.6 Mq of small clus- 
ters, while its 20 Mq counterpart produces ~ 0.103 Mq. Our results point to smaller 
amounts of dust formed in the ejecta of Pop. Ill supernovae by a factor ~ 5 compared 
to values derived by previous studies, and to different dust chemical composition. Such 
deviations result from some erroneous assumptions made, the inappropriate use of clas- 
sical nucleation theory to model dust formation, and the omission of the synthethis of 
molecules in supernova ejecta. We also find that the unmixed ejecta of massive Pop. Ill 
supernovae chiefly form silica and/or silicates, and pure silicon grains whereas their 
lower mass counterparts form a dust mixture dominated by silica and/or silicates, pure 
silicon and iron sulphides. Amorphous carbon can only condense in ejecta where the 
carbon-rich zone is deprived of He + via the nucleation of carbon chains and rings char- 
acteristic of the synthesis of fullerenes. The first dust enrichment to the primordial gas 
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in the early universe from Pop. Ill massive supernova comprises primarily pure silicon, 
silica and silicates. If carbon dust is present at redshift z > 6, alternative dust sources 
must be considered. 

Subject headings: astrochemistry — supernovae: general — early universe — molecular 
processes — dust 



1. Introduction 



For the last two decades, many efforts have been developed to try and understand the role of su- 
pernovae (hereafter, SNe) as potential dust contributors to the Universe. On the observational front, 
the first evidence for dust formation in a SN event was observed during the explosion of the Type 
Hp supernova SN SN1987A. Its extensive observational coverage at mid-infrared (IR) wavelengths 
revealed the presence of the fundamental and overtone transitions of a few molecules, specifically 
CO and SiO, as early as ~ 120 days post-explosion, and the formation of dust grains after day 



400 Meikle et al. 



1989 



Moselev et al.lll989l ; iRoche et al 



1991 



Danziger et al. 



1991 



Wooden et al 



19931 ). The exact composition of these condensates is yet not known exactly but recent modeling of 
the IR excess mea sured by Spitzer usin g a clumpy ejecta proposes a mixture of amorphous carbon 
and silicate dust (jErcolano et al.l 120071 ) . Several other galactic and extra-galactic SN explosions 
have bee n monitored since t hen, and th e detection of CO, SiO and dust in their ejecta has proved 
positive (IKotak et al J 12005 . 120061 . 120091 : ISugerman et al.1 120061 ) . However, wether a SN explosion 
produces large amounts of dust or not in its ejecta is still a matter of debate. Indeed, the values of 
the dust mass synthesized derived from observational IR data are never larger than ~ 10~ 2 Mq. 
Large amounts of dust (> 10 8 Mq) at high redshift (z > 6) have been conjectured to explain the 
sub-milli meter observations of distant quasar s and metal measure ments in Damped-Lyman- Alpha 



Systems (jPei et al.l Il99ll ; iPettini et al. 



19941 ; iBertoldi et al.l |2003j) . In view of the small amount 



of dust detected in local SN ejecta, these findings pose problem for a dust SN origin in the early 
universe and the po ssible synthesis of solids in the explosive ends of massive Population III stars 
(hereafter Pop. Ill, (|Dwek et al.ll2007r A. 



The formation of dust in SN events has also been tackled theoretically and several studies 
have been conducted for local and high redshifted SNe. Two approaches are usually used for 
modeling purposes: (1) the thermodynamical equilibrium (TE) assumption, and (2) the classical 
nucleation theory, or CNT. The first approach is assumed in models aiming at deriving dust and 
mixing properties in local SN ejecta from the study of isotopic anomalies in meteorites (jLodders 



1996 



Travaglio et al.lll999h . The assumption of TE implies chemical equilibrium. The derivation 



of condensation temperatures and chemical compositions in a pressure-temperature equilibrium 
phase diagram for the various solids under study may be adeq uate to describe the formation of 
molecules in stellar phot ospheres where high pressures prevail (ITsujiJ Il973l : ICherchneff &; Barker 



1992 



Helling et al.lll996l ) but is clearly inappropriate for lower density circumstellar environments, 
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in particular the dynamical flows of SN ejecta. As we will see below, the chemistry is neither at 
equilibrium nor at steady state over the timespan under which dust condenses in SN ejecta, while 



used bv various groups for low and high redshift SNe (Kozasa et al. 


1989 
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Clavton et al. 


1999: 


Todini & Ferrara 


2001 


Nozawa et al. 


2003: 


Schneider et al. 


2004 


Bianchi & Schneider 


2007 


). CNT 



was first developed to explain the forma tion of water droplets in the Earth atmosphere under 



equilibrium conditions (IFeder et al. 



1966) and used i n initial attempts to model dust formation 



in the interstellar medium (jDraind Il979l ). However, iDonn &: Nuthl (|l985l ) first objected to the 
appropriateness of CNT to describe dust synthesis in space, i.e., when the formation of solids from 
the gas phase occurs in low-density systems far from equilibrium, e.g., in circumstellar environments 
such as stellar winds of evolved stars or SN ejecta. Indeed, CNT involves the calculation of critical 
cluster sizes for nucleation and bulk material concepts like specific surface energies used in the 
derivation of the free energy barrier for nucleation. However, the derived critical clusters are 
usually on the atomic scale, making clear that CNT bulk property concepts are arguable for atoms 
and small molecules. The existence of a steady-state critical cluster distribution is questionable 
if the molecular phase from which those clusters form is not at steady state. Finally, even in 
high-pressure experiments where equilibrium conditions hold, the nucleation steps involved and 
the nucleation products are different from those predicted by CNT. These limitations have been 
recognized for some time by the various scientific communities modeling the synthesis of ceramics 
or soot in flames in the laboratory. In such experiments, the formation of dust proceeds through 
a two-step mechanism involving the nucleation of small molecular clus t ers and the growth of those 
clust ers from surface deposition and coagulation (jMcMillin et al.lll996l ; |Pratsinislll998l : IWooldridge 
19981 ). All those processes are kinetic in essence and usually occur far from equilibrium. 



A stochastic, kinetically-driven approach must thus be used to describe dust formation in space, 
in which gas-phase molecules transform into small dust clusters via chemical kinetic processes with 
no prior assumption on the nature of solids that should condense from the initial gas mixture 
under study. A major drawback is the partial lack of information on the different steps involved 
in nucleation and their related chemical reaction rates, and on the intermediate species formed 
during coagulation. However, many recent studies of ceramic synthesis in flames highlight the 
role of key species and provide new insights on these matters. A first attempt to model the 
formation of dust in SN ejecta at high redshift using a chemical kinetic approach based on a careful 
descrip tion of the formation of m olecules and the nucleation of small solid clusters was carried 
out bv ICherchneff &: Lilly! (|2008l ). In a more complete study with revised chemical processes, 
Cherchneff &: Dwekl ()2009l . hereafter CD09) show that Pop. Ill SNe are efficient at forming key 
molecules such as O2, CO, SiS and SO in substantial amounts (up to ~ 35 % of the progenitor 
mass). Those molecules have a crucial impact on the dust synthesis as a) they deplete the gas 
phase from potential heavy elements that could otherwise be included in grains and b) they act as 
coolants to the local gas through emission in their ro-vibrational transitions. For example, O2 and 
CO are the prevalent molec ular coolants through their rotational transitions during the collapse 
phase of molecular clouds ([Goldsmith &; Langerl Il978l ) . CD09 also show that the formation of 
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chemical species in the ejecta is neither at equilibrium, nor at steady state, a common assumption 
used in all existing studies. 

In the present paper, we study the formation of small molecular clusters involved in the nu- 
cleation step of the synthesis of dust in the explosive ends of Pop III stars. Those clusters are of 
paramount importance in the determination of the dust chemical composition as well as the final 
dust mass formed as they are the bottleneck to condensation processes. We consider small clusters 
involved in the condensation of several solid compounds commonly found in evolved, circumstellar 
environments and whose chemical compositions are in agreement with the elemental composition 
pertaining to SNe. They include pure metallic grains, metal oxides, metal sulphides, silica and 
silicates, silicon carbides and amorphous carbon dust. We study their formation for the two mixing 
cases described in CD09, i.e., fully microscopically- mixed and unmixed ejecta, for a zero-metallicity, 
very massive progenitor exploding as a pair-instability SN (or PISN) and a low-mass progenitor 
exploding as a core-collapse SN (CCSN). In §2, we fully describe the chemistry of cluster formation 
from the gas phase. §3 briefly presents the SN models assumed in this study whereas results are 
presented in §4 along with a critical assessment of existing studies, and an analysis of key parame- 
ters to the dust formation processes. Finally, a summary of the results and their interpretation in 
the context of the early universe are presented in §5. 



2. THE CHEMISTRY OF DUST NUCLEATION 

Dust formation in circumstellar environments does not occur under equilibrium or steady 
state conditions due to the complexity of the dynamics and the physics involved in such media. As 
in the laboratory, dust condensation and growth occur through stochastic, kinetically-controlled 
processes. For example, in Asymptotic Giant Branch stars, the chemistry characterizing the gas 
layers located above the stellar photosphere and where dust forms is not at equilibrium despite the 
large gas densities an d temperatures en countered, owing to the action of periodic shocks induced by 



the stellar pulsations (|CherchnefHl2006l ) . In supernova ejecta, the passage of the blast wave through 
the helium core and the deposition of radioactive energy in the ejecta trigger an active chemistry 
far from equilibrium or steady state, as shown in CD09. As such, the formation of dust resembles 
that occurring in pyrolysis, flame, or vaporization experiments, and consist of two phases: (1) the 
nucleation of small molecular clusters from the gas phase determined by chemical kinetic processes; 
(2) the growth of those clusters to dust grains from stochastic coagulation and surface addition 
processes. 

Depending on the initial composition of the gas, the clusters are of different chemical nature. In 
an oxygen-rich environment, metal oxides and various silicates are expected to form whereas a large 
carbon content fosters the formation of metal carbides and various forms of solid carbons (polycyclic 
aromatic hydrocarbons, i.e., PAHs, fullerenes, carbon chains and rings etc.). Clusters of pure iron, 
silicon, and magnesium also form in vaporization experiments as well as metal sulphides, when 
atomic sulphur is present. We will describe below the various clusters considered in the present 
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study and their related formation chemistry. Those clusters are the building blocks of the various 
'astronomical' types of dust observed in circumstellar environments and are all synthesized in flame, 
condensation, and vaporization experiments in the laboratory. The nucleation stage described in 
this study involves the formation of ID- molecular chains, 2D- rings and 3D- small cages up to a 
maximum of 10 atoms. This cluster size upper limit is somehow arbitrary and is motivated by the 
impact of the cluster geometry on the growth process, i.e., a cluster with a 3D structure like a cage 
or a twisted ring can be seen as a small seed nucleus for coagulation with a larger collision cross 
section than those of chains or rings. The clusters and the dust that condenses out of them are 
listed in Table 1. 

All chemical pathways leading to the formation of chains include neutral-neutral processes like 
bimolecular and radiative association reactions whereas destruction is described by thermal frag- 
mentation, reaction with Compton electrons, and neutral-neutral processes (i.e., oxidation reactions 
for carbon chains). In the special case of rings, their structure precludes a large chemical reactivity 
and we assume that the polymerization of the gas-phase monomer proceeds via neutral-neutral pro- 
cesses while the destruction of polymers results from thermal fragmentation only. Hence, reactions 
of atomic oxygen and He + ions with rings are ignored. All processes are discussed in detail below 
for the various clusters and a list of all chemical reactions involved in the nucleation of clusters 
and their respective rates is given in Table 2. The full chemical network considered in this study 
includes Table 2 and Table 5 of CD09, the latter describing the chemistry involved in the synthesis 
of molecules in SN ejecta. As all those chemical processes apply to hot (~ 2000 K) and dense 
(~ 10 9 mol. cm 3 ) gas, these two tables represent an extremely useful compilation of reactions and 
rates for modeling the dust chemistry of other circumstellar environments such as Type Ic super- 
novae, Novae, the shocked inner regions of AGB stellar envelopes, Supergiants and LBVs outflows, 
the wind of RCrB stars, or the colliding winds of Wolf-Rayet stars. 



2.1. Oxygen-Rich Clusters 



Amorphous and crystalline silicates are the most abundant family of solids found in space. 
Olivine grains including forsterite (Mg2Si04) and fayalite (Fe2Si04) are observed in the wind of 
oxygen-rich low-mass stars ascending the A GB when pyroxene dust (enstatite MgSi03 and ferros- 
tite FeSIOs) are found around young stars (IMolster fc Waters! 120031 ) . For silicates, no gas phase 
monomer precursors exist from which solids can grow by direct coagulation. Despite the wealth 
of studies dedicated to the formation of silicates and ceramics in the laboratory, the identifica- 
tion of molecular precursors to silicate formation is still a difficult task. Flame aerosol technology 
experiments high light the important role of heterogeneous seed cluster s in the formation of sili- 



cate compounds (jMcMillin et al 



1996 



Pratsinis 



1998 



Fe/Mg/Si/02-rich vapor mixture in condensation experiments, (jRietmeiier et al 



Wooldridgelll998l). In presence of an initial 



1999 



2002J) show 



that the final products include eutectic Mg-SiO magnesiosilica or Fe-SiO ferrosilica condensates but 
no mixed ferromagnesiosilica grains. The formation of magnesio- and ferrosilica involves the coales- 
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cence of heterogenous small clusters of silica, magnesium oxides and iron oxides, and pure metallic 



clust ers dependent of the initial ratio of metals to gas-phase SiO (jKaito et al.ll200a : iKamitsuii et al 



20051 ) . We discuss below the possible reaction channels to the formation of small silica an d metallic 



oxide clusters. Clus ter structure and stability have been derived for silicon oxide clusters (ILu et al 



20031 ). metal oxides (IZiemann et al.lll9 90; Kh oler et aUll997l : lBhatt & Fordll2007l : iReillv et alJl2007l ). 
and pure metallic clusters (jjeilinek k, Aciolill2002l ). The structures for (SiO)„, (MgO) m and (Mg) m 
are schematically represented in Figure 1 for n = 1 — 5 and m = 1 — 4, respectively. 



2.1.1. Silica 



In space, the 20 ^m band of Si02 was detected in Herbig Ae/Be stars with ISO (jBouwman et al 



2001 



) , while the mode ling of Cas A IR spectra from Spitzer suggests Si02 as a potential dust candi- 



date (|Rho et al.ll2008l ). In the laboratory, silica clusters form following two r outes depending on the 



hydro gen content of the gas from which they nucleate. Where H is present, IZachariah Burguess 



(|1994l ) report the importance of the OH hydroxyl radical in the formation of silica clusters in 
H/O/Ar-rich flames fueled with silane (S1H4), using Laser-Induce Fluorescence techniques to char- 
acterize the gas phase species entering in the conversion of gas to particle processes. Indeed, OH 
triggers the formation of gaseous Si02 through the reactions (see CD09 for more detail) 

SiO + OH — ► Si0 2 + H, (1) 

SiO + H 2 — > Si0 2 + H 2 , (2) 
and consequent cluster growth occurs from the Si02 polymerization reaction 

Si0 2 + (Si0 2 ) n — > (Si0 2 ) n+l . (3) 

where 1 < n < 4. 

When the gas is hydrogen-poor, the nucleation of silica occurs through direct polymeriza- 
tion of SiO to form small (SiO) n molecular clusters, precursors to solid silicon monoxide SiO. At 
temperatures > 900 K, solid SiO is unstable and directly decomposes to silica Si02 and solid Si. 
We thus consider that half of the total number of (SiO)s particles will be in the form of (SiOa)5, 
whereas the other half will be in solid (Si)s- It is not clear wether this disproportionation effect 
gives rise to two separated condensate populations o f SiO? and pure Si o r to solids made of pure 



Si clusters embedded in an amorphous silica matrix (IMamiya et al.l 120011 ). We neglect this aspect 
in the present study and consider that solid SiO gives rise to the formation of two types of solids, 
whatever their final form. SiO is first synthesized from the Si oxidation reaction (see CD09 for 
more detail) 

Si + 2 — > SiO + O, (4) 

by reaction with CO 

Si + CO — > SiO + C, (5) 
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and the radiative association reaction 



Si + O — > SiO + hu, 



followed by the polymerization of SiO 



SiO + (SiO) n — > (^o) n+1 . 



(6) 



(7) 



with 1 < n < 4. In pyrolysis experiments, reaction [7] usually happens at almost collision rate and is 
at or near high-pressure limit at room temperature. However, as the temperature increases, the rate 
decreases due t o the onset of energy transf er effects. The rate is also highly pressure-dependent and 
is estimated bv lZachariah &; Tsana (|1993l . hereafter ZT93) for the high temperatures encountered 
in the combustion flames and applicable to SN ejecta. For our ejecta, we scale the rates derived 
by ZT93 for 1 atm pressure according to the total number of collisions per cm -3 for the ejecta 
pressure under consideration. The resulting rates for the polymerization reactions of type reaction 
[7] are a factor ~ 100 smaller than those of ZT93 at 1 atm pressure. We model both the formation 
of (SiC>2)n and (SiO) m up to n = 4 and m = 5. The most stable form of (SiC>2)4 has a rhombus 
chain structure with double oxygen bridged whereas the most stable s tructure of (SiO )4 is a buckled 
eight-membered ring and that of (SiO)s has a double ring structure (|Lu et al.ll2003l ). 



The destruction process for those silica clusters is the thermal fragmentation reaction 



(SiO) n + M — > (SiO) m + (SIO) n _ m + M, (8) 

where M represents the ambient gas , and 2 < m < 5. Theoretical estimation of dissociation energies 
for silica clusters bv lLu et al.l (|2003l ) shows that the most common fragmentation products are gas- 
phase SiO, (SiO)2 and (SiO)3, as found in experiments on silica cluster material. This result implies 
that those species must be rather stable in the gas phase, and for the sake of simplicity, we assume 
a rate for reaction [8] equals to the rate for the thermal fragmentation of SiO (see Table 9 in CD09). 



2.1.2. Iron Oxide 



Iron oxide has bee n proposed as a potential contributor to the 23 iim band s in young stars 



Rho et al 



Bouwm an et al.ll2000l ). in the binary AFGL 4106 system (jMolster et al.1 119 99). and in Cas A 



2008). The reaction of iron with oxygen to generate iron oxides is at the basis of funda- 



mental processes in environmental corrosion, biological oxygen transport or oxide film formation. 
Of special interest to the present study is the gas phase iron monoxide formation reaction 



Fe + 2 — ► FeO + O, 



(9) 



whose rate was measured at high temperatures bv lFontjin et al.l (jl973l ) and lAkhmadov et al.l (|1988l ). 
The following reaction with carbon dioxide also occurs 



Fe + C0 2 — ► FeO + CO, 



(10) 
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and its rate was measured by ISmirnovl (|2008l ). Finally, we consider destruction of iron monoxide 
through thermal fragmentation and its reaction with atomic oxygen 



FeO + O — > Fe + 2 , 



(11) 



for which a rate was derived by ISelf fe Pland (j2003l ) in their study of iron oxide formation in 
the Earth upper mesosphere. The thermodynamically stable neutral clusters resulting from laser 
ablat ion of an iron rod bathed in an oxidizing gas are of the form (FeO) m for m < 10 (jShin et al. 
2004). We thus define the nucleation mechanism to small (FeO) m clusters in the gas phase by 
successive addition of the FeO monomer 



FeO + (FeO) r 



(FeO) n 



(12) 



with 1 < n < 3. The rate for this reaction is chosen similar to that of reaction [7] Destruction of 
small iron oxide clusters proceeds through thermal fragmentation such as reaction Iron oxide 



(Wang et al. 


1996: 


Reillv et al. 


2007) 



2.1.3. Magnesium oxide 



Magnesium oxide MgO (periclase) was proposed as a component of the dust around the low- 
redshift quasar PG 2112+059 for its contribution to the 20 fim spectrum region (jMarkwi ck-Ke mper et al 



2007) and was suggested as potential seeds to silicate condensation in O-rich AGBs ([Bhatt fe Ford 



20071 ). The nucleation of magnesium oxides involves the formation of MgO molecules in the first 
place through reaction 



and 



Mg + Q 2 — > MgO + O, 



Mg + C0 2 — ► MgO + CO, 



and MgO destruction by thermal fragmentation and reaction with atomic oxygen 

MgO + O -^Mg + 2 . 



(13) 



(14) 



(15) 



The rate for reaction [13] was estimated by iKhashireninov et al .1 (j 19821 ) whereas those of reactions 
[T4l and [15] were assumed equal to the equivalent reactions involving atomic Fe. We follow the 
formation of small magnesium oxide clusters (MgO) n with n = 1 — 4. As seen from Figure 1, the 
{MgO) 2 and (MgO)s clusters have a rhomb us and triangular structures, respective ly, whereas the 
lowest energy structure of (MgO)^ is cubic (jKholer et al.lll997l ; iBhatt fe Fordll2007l ). We therefore 
describe MgO cluster formation by the direct polymerization of MgO according to reaction 



MgO + {MgO) n — ► {MgO) n 



+1) 



(16) 



where a rate similar to that of reaction [7] has been assumed. Small periclase clusters are destroyed 
according to the thermal fragmentation reactions similar to reaction [8] 
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2.1. 4- Alumina 



Alumina AI2O3 has been detecte d through its 11 fim band in evolved, oxygen-rich AGB stars 



( Stencel et al. 



1990 



Cami et al.ll2002l ) and was fo und in pre-solar me teorites with an isotopic com 



position pointing to its O-rich AGB stella r origin (INittler et al.lll997h. M odels of the mid-IR excess 



of the low-redshif t quasar PG 2112+059 (jMarkwick-Kemper et al 



2007), the Blue Luminous Vari- 



able Eta Carinae (jChesneau et al.ll2005l ). and the SNR Cas A (IRho et al.ll2008l ) also use alumina as 
a potential dust component. In the laborarory, laser-induced alu minum evaporatio n in an helium 
bath seeded with oxygen produces small aluminum oxide clusters (IDesai et al.lll996l ). Small cluster 
molecules like A^O^ (x = 3 — 5) show a rhombus structure, including AI2O2 rings. The chemical 
pathway to those small clus ters involves the fo rmation of the molecule AlO, and its further reac- 
tion with oxygen and itself (jCatoire et al.ll2003l ). In our SN ejecta, we only consider the formation 
of AlO as aluminum has an initial low mass abundance compared to other metals reacting with 
oxygen (see Tables 2 and 3 of CD09). Furthermore, the formation routes to alumina are still very 
speculative and we prefer deriving AlO mass yields as a first indication of the alumina content of 
the ejecta. AlO synthesis occurs through the following reactions 



Al + 2 — ► AlO + O, 

Al + O + M — ► AlO + M, 
Al + C0 2 — ► AlO + CO, 
and finally, in the presence of hydrogen, 

Al + H 2 — ► AlO + OH. 



(17) 

(18) 
(19) 

(20) 



Reaction rates for reactions 



({2003] ), and iMcClean et al 



ana 

19931 



121031 CHI and[20]were derived bv lGarland & Nelson! (119921 ) . iGatoire et al 



(119931 ) . Destruction process for gas-phase AlO molecules are thermal 



fragmentation by the ambient gas, collisions with Compton electron, and reaction with He + . 



2.2. Metal sulphides 



The metal sulphides MgS and FeS have been proposed as carriers for several IR bands observed 
in evolved, carbon-rich, circumstellar environments. Specifically, MgS is most likely responsible for 
the broad 30 Atm emiss i on band observed in carbon stars, proto-planetary and planetary nebulae 
(jBegemann et al.1 11994 : IVolk et al.l 120 02). while FeS has been proposed as the carrier of the 23 
/im feature observed in two planetary nebulae with ISO (jHony et al.ll2003l ). On Earth, the iron 
sulphide clusters (FeS)2 and (FeS)4 act as evolutionary ancient prostheti c groups in sustaini ng 
fundamental life processes such as biological electron transport in proteins (j Johnson et al.ll2005l ). 



As no reaction between metals and sulphur is documented, we use the S and O isovalence to 
assess reaction rates for metal sulphide formation. Following reaction \13\ MgS formation starts 
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with reactions 



and 



Mg + S 2 — > MgS + S, 



Mg + SO — ► MgS + O. 



(21) 



(22) 



to both of which we assign the rate of reaction [13] Destruction of MgS occurs via thermal frag- 
mentation and the reverse processes of reactions [21] and [22] with assigned rates similar to those for 
reaction [15] Finally the clustering of MgS proceeds via the reaction 



MgS + (MgS) n — > (MgS) n 



+1, 



(23) 



where a rate similar to that of reaction [7] has been assumed. The ground state structures of MgS 
clusters are as follows: a rhombus structure for (MgS)2, a planar triangular structure for (MgS)3, 
and a distorted cubic structure for (MgS)4. 

Similar reactions and rates are considered for the growth of troilite FeS clusters. Formation 
processes are 

Fe + S 2 — > FeS + S, (24) 



and 



followed by 



Fe + SO 



FeS + O, 



FeS + (FeS) n — ► (FeS) n+1 . 



(25) 



(26) 



Destruction of FeS occurs via thermal fragmentation and the reverse processes of reactions [24] and 
[25l Clusters of FeS have the following structures: rhombic for (FeS)2, planar and triangular for 
(FeS)3, and distorted cubic for (FeS)4. 

We do not consider the formation of solid SiS in the present study, although gaseous SiS being 
isovalent to SiO, similar nucleation trends to (SiO) n synthesis could be expected. Solid SiS is known 
to disproportionate into solid SiS2 and solid Si upon annealing at temperatures greater than 900 
K (IByerley fc Teolll973l ). SiS2 solid was propo sed to explain the 21 //m emiss i on ba nd observed in 
several carbon-rich Proto-Planetar y Nebulae (Goebe] 1993 : Begemann et al. 1996 ). although the 
identification is still controversial ( Zhang et al. 20091 ). However, the large abundances of gaseous 
SiS predicted to form in SNe ejecta by CD09 could spawn the synthesis of solid silicon disulphide 
in SNe. 



2.3. Pure Metal Clusters 



There exists no observational evidence for the existence of pure metal clusters in evolved 
circumstellar environments. However, iron clusters have been proposed as a comp onent of dust 
in O-rich AGB stars to explain the required near-IR opacity ([Kemper et al.l 12002 ) . They were 



-11- 



also included in dust models aiming at reproducing the IR excess of Cas A observed by Spitzer 
(|Rho et al.ll2008l ). Pure iron, magnesium and silicon clusters are easily formed in the laboratory. In 
particular, Fe nanoparticles gas-phase synthesis from iron pentacarbonyl Fe(CO)s decomposition 
has been e xtensively studied as iron particles are important catalysts in the formation of carbon 
nanotubes (jWen et al.l 120071 ). The Fe small cluster nucleation proceeds according to the following 
steps 

Fe + (Fe) n + M — > (Fe) n+1 + M, (27) 
where n = 1,3. The coagulation of small clusters is also considered 



Fe 2 + Fe 2 



Fe A 

Fe 3 + Fe, 



(28) 



when destruction occurs through the thermal fragmentation reaction 



Fe n + M — > Fen-! + M. 



(29) 



The rates for reactions 1271 [28l and 1291 have all been estimated by (jGiesen et al.ll2003l ). Similar 
processes and rates have been assumed for the nucleation of pure Si and Mg clusters. 



2.4. Carbon-Rich Clusters 

Along with silicates, carbon-rich clusters like amorphous carbon (AC) and silicon carbide (SiC) 
are ubiquitous condensates in space. In the present study, we consider the nucleation of both AC 
and SiC clusters in primordial SN ejecta. 



2.4-1- Amorphous Carbon 



The main providers of AC grains in galaxies are carbon-rich AGB stars in the late stages of 
their evolution, but other formation loci include RCrB stars, the colliding-wind region of carbon-rich 
Wolf-Rayet binary systems, and SN ejecta. Whenever the carbon content of a gaseous medium is 
larger than that of oxygen, amorphous carbon dust may condense. The nucleation to solid clusters 
i nvolves different carbon routes dependent of the temperature and the hydrogen content of the gas 
(|cherchnenil2009l ). When hydrogen is present and for the low temperatures encountered in premixed 
acetylene flames (T gas ~ 1500 K), unsaturated hydrocarbons form, among which acetylene C2H2 
and its isomeric radical vinylidene (HHCC), the ethynyl radical C2H, the propargyl radical C3H3, 
and vinylacetylene Ciji\. The formation of the planar, aromatic ring benzene CqHq results from 
successive additions of acetylene and its isomer, or the direct reaction of two propargyl radicals 
([Miller &: Melliud 120021 ) . Its further growth to Polycyclic Aromatic Hydrocarbons (PAHs) occurs 
through successive addition of acetylene on a romatic radical sites (jFrenklach &: Feigelsonl [1989J; 



Cherchneff et al. 



1992 



Richter fc Howard! 12002 ). 
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When the environment is hot {T gas ~ 4000 K) and hydrogen- free, the nucleation of carbon 
clusters involves the synthesis of pure carbon chains of sp hybridization up to Cg. For Cio, the 
most stable structure is ringed and the closure of the chain to large monocyclic rings occurs. 
The coalescence of such rings in bicyclic and tricy clic rings as precursors to fullerene formati on is 
observed in pulsed laser vaporization experiments (lyon Helden et al. 119931 ; iHunter et al.lll993l ) and 
the condensate is in the form of fullerene-like soot ( Jaeer et alJboOfll ). 



For the purpose of this study, we explore the kinetics of the formation of carbon chains up 
to Cg and the closure to the monocyclic ring C\q. As our ejecta are considered hydrogen- free as 
hydrogen is not microscopically mixed with other constituents, we disregard the formation routes to 
aromatics involving hydrocarbon s. Our kinetic network for chain formation inclu des the radiative 



association reactions of the type (jClayton et al. 



1999 



2001 



Deneault et al.ll200fil ). 



C + C n ^ C n+ i + fw, 



(30) 



with n = 1 - 
(f m =lx 10- 17 
(ten] = 1 x io- 10 



- 9. 

cm 3 s~ 
cm 3 s~ 



1, 



Andreazza fc Singhl (j 19971 ) derived a low associative rate 



For n 

1 ) whereas the rate for larger chain synthesis is assumed much faster 



as larger chains are able to qu ickly stabilize their formation com- 
plex owing to their numerous vibrational degrees of freedom (jClayton et al.lll999l ). An additional 
gr owth process for c hains is their mutual collision and binding by edge carbon atoms. According 



to ISchweigert et al.l ()1995l ). such mechanisms have energy barriers less than 0.5 eV. The thermody- 
namics of the addition react i ons ar e estimated using thermodynamic data for small carbon chains 
calculated by IClavton et al.l (|200ll ) and we find that these additions are all exothermic processes. 
Therefore, we include the additional growth pathways to chain synthesis 



Cn—m + C m y C n , (31) 

with m = 2 — 7. The rate for reaction [3D is assumed to be collisional with an estimated value of 
1 x 10~ 10 cm 3 s -1 . 

The dominant destruction channels of carbon clusters included in the model are firstly, the 
thermal fragmentation of the chains through collision with the ambient gas 



C n + M — > C n - m + C m + M, 



(32) 



secondly, the oxidation of the chains through reaction of atomic oxygen with end-cap carbons to 
form CO 

O + C n — ► C n -i + CO, (33) 
thirdly, the collision with Compton electrons 



C n + e Compton 

and fourthly, the destruction by He + 

C n + He^ 



Cn—m C m ~\~ 6 



Compton ' 



C n -! + C+ + He. 



(34) 
(35) 
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Th e rates for rea c tions [32] are those used for C2 in CD09, while the rate for reactions [33] is 
that of Clayton et al. ( 1999 . 12001 ). and finally, the rate for reactions [M] is derived following the 
formalism explained in CD09. 



2.4.2. Silicon Carbide 



Silicon carbide grains are ubiquitous in various evolved circumstellar environments (IMolster fc Waters 



2003). In AGB stars, the SiC 11.3 um band is observ ed either in emission or absorption for low and 
high mass loss rates, respectively ([Speck et al.lll997l ). but there is no detection of this band in SN 
ejecta to date. Presolar SiC grains bearing the isotopic signatures pertaining to SNe are however 



identified in meteorites (Bernatowicz et al 



19871 ). In SN ejecta, the formation of SiC grains implies 



the mixing of Si- and C-rich non-adjacent layers, a situation not found in the dust formation layers 
of car bon stars where SiO and C-b earing chemical species are simultaneously present in large quan- 
tities ([Willacy &i Cherchnefilll998l ). In the laboratory, SiC whiskers and nanorods are produced by 
different methods from vaporization of SiC bulk material, to the surface reaction of gas-phase SiO 
on carbon fil ms or nanoclusters, or th e gas phase formation from chlorosilanes, carbon monoxide, 
and methane ([Zhou &: Seraphinlll994[ ). In the present study, we consider the formation of SiC and 
(SiC)2 according to the following reactions: 



and 



followed by the nucleation reaction 



Si + CO — > SiC + O, 
Si + C 2 ^ SiC + C, 
C + SiO — >• SiC + O, 

C + Si — ► SiC + hv, 
SiC + SiC — > (SiC) 2 . 



(36) 
(37) 
(38) 

(39) 
(40) 



The rate for reaction 00] is assumed to procee d with the rate of reaction [7] The (SiC)2 cluster 
has a favored rhombic structure (jYadaf et ajj[2006j ) similar to that of (SiO)2 in Figure 1, but a linear 
structure is also possible. The destruction channels for SiC are reactions with He + , the reverse 
processes of reactions I3"6ll38[ and thermal fragmentation, while the latter is the only destruction 
channel considered for (SiC)2- 



2.5. From Clusters to Dust: Considering Fe and Mg depletion 



The condensation stage in dust synthesis involves the coagulation of clusters and the growth 
of seeds through additional adsorption and chemical reactions on the grain surface. Those pro- 
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cesses are extremely complex and their description is out of the scope of this paper. However, 
we can rely on useful information coming from the synthesis of smokes in the laboratory to make 
prescriptions for the dust formed from the various types of clusters synthesized in SN ejecta. We 
will adopt the following approach for the derivation of the dust upper limits in the rest of this 
paper. We may encountered situations where large amounts of (SiO)s rings are synthesized in a 
medium rich in atomic Mg and/or Fe. These SiO rings will lead to the formation of metastable 
solid silicon monoxide which disproportionates into a mixed phase of pu re Si nano-crystals embed 
ded in amorphous s ilica clusters at temperatures larger than ~ 900 K (IFiiglein fe Schubert 



1999 



Kaplatis et al.ll2005l ). According to §2.1, such a rich Mg- or Fe-SiC^/Si mixture does not lead 
to the formation of ferroma gnesiosilica, but to separated populations of metastable ferrosilica o r 



magnesiosilica compounds (jRietmeiier et al.l 119991 . 120021 ; iKaito et al.ll2003l ; iKamitsuii et al.ll2005l ). 
Furthermore, the formation of additional magnesium silicide Mg2Si and periclase MgO clusters is 
observed in the Mg-rich zone of the smoke resulting from the combined evaporation of Mg and SiO 



solid s whereas pure Si and f orsterite Mg2SiC>4 clusters are observed in the Si-rich zone (|Kaito et al 
200.4 iKimura fc Nuthlbood ). 



In view of the complexity of the products and because the exact composition of those metastable 
compounds is not truly characterized, we here consider two simple cases, once the nucleation of 
molecular clusters has been completed: (1) a case for which all pure Fe and/or Mg atoms and clus- 
ters are depleted in the dust condensation process leading to either Mg2Si04 (forsterite) or Fe2SiC>4 
(fayalite) dust; (2) a case for which no Fe/Mg depletion occurs. For the sake of simplicity, we ignore 
the possible formation of magnesium silicide Mg2Si and pyroxene compounds like enstatite MgSi03 
and ferrosilite FeSi03. The former case is somehow extreme assuming a depletion efficiency during 
the grain formation process of 1 for the metals Fe and/or Mg, whereas the latter may depict a 
situation more in tune with the ejecta of local SNe. Indeed, in the CCSN SN1987A, optical lines of 
[Mgl] and [OI] were observed to decline with time according to the ejecta optical depth whereas the 
[Sil] line was observed to have a sharp d ecline owing t o Si de pletion at the onset of dust formation 
(jLucv et al.l 1 19891 ; IWooden et al.l Il993l ) . iMeikle et al.l (| 19931 ) also observed a sharp decline in the 
1.2 fim and 1.645 /jm [Sil] lines at the onset of dust formation compared to the 1.257 /im [Fell] 
line. Observations thus indicate that Mg and Fe are not as depleted as Si, and it should be borne 
in mind that reality probably lies within these two extreme cases for metal depletion. We then 
derive upper limits for Mg2Si04 and Fe2Si04 according to these two cases. To satisfy stochiometry 
and form forsterite and fayalite and because of the large dioxygen ejecta content at late times af- 
ter explosion, we assume that the magnesiosilica and ferrosilica pre-seeds will be oxidized through 
O2 adsorp tion and decomposition on the grain surfaces at the low temperatures encountered at 
day 1000 (jOver fe Seitsonenl l2009l ) . This is a plausible mechanism for oxidation as O2 forms at 
t > 700 days after the decrease of He + abundance as seen in CD09 and §4. The oxidation process 
will partly deplete the O2 reservoir, resulting in smaller O2 mass yields ejected at day 1000 when 
metal depletion is considered. In case no metal depletion and formation of forsterite or fayalite are 
assumed, the major SiO-based solid in the ejecta is silica Si02- 
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THE EJECTA MODELS 



For the present study, we pay attention to two Pop. Ill progenitor masses, a 170 Mq PISN, 
and a 20 Mq CCSN. We consider two mixing conditions: (1) a fully mixed ejecta in which all 
elements have been microscopically-mixed in the helium core; (2) an unmixed, stratified ejecta in 
which each layer elemental composition of the helium core reflects the prior nucleosynthesis stages 
of the progenitor, except for the inner most mass zone whose composition is the result of explosive 
nucleosynthesis and radioactive decay. None of these scenarios is fully satisfactory in term of ejecta 
description, but a stratified ejecta made of fully microscopically-mixed zones of various chemical 
compositions seems more appropriate to describe the chemistry of the gas filaments and clumps 
resulting from the explosion. However, we consider fully microscopically-mixed ejecta as those 
are considered in most of the existing studies. Mixing is a key parameter to the ejecta chemistry 
and we will discuss it s impact on dust formation in §5, in the light of 2D/3D explosion models 
(|Frver Warrenll2002l ; iKifonidis et al.ll2003l . 120061 ; lHammer et alj|2009l ). A complete description of 
the ejecta models is given in CD09, in which full details on the physics and chemistry at play and 
the model parameters used in the present study can be found. 

As a brief su mmary, the g a s par ameters like temperature and number density are derived 
from the study of iNozawa et al.l (|2003l . hereafter NK03) for both mixed and unmixed ejecta, and 
include the radioactive decay of 56 Ni, 56 Co, and 56 Fe which induces the creation of a population 
of Compton electrons in the gas through 7-rays down-scattering. The impact of the ultra violet 
(UV) radiation field resulting from this 7-rays degradation on molecules and dust clusters was 
assessed by CD09 who found that the destruction of molecules and dust precursors by this UV 
radiation field was not important. We thus ignore UV radiation for the rest of the present study. 
The c hemical compositions for fully microscopically-mixed ejecta are those of lUmeda &: Nomoto 
(J2002J, hereafter UN02) when elemental mass yields for unmixed ejecta are taken from NK03. The 
chemical compositions of our present models and the zoning considered in the unmixed ejecta are 
given in Tables 4 and 5 of CD09. 

The gas parameters of fully mixed ejecta over the timespan of interest (100 to 1000 days) are 
summarized in Table 3 for the two progenitors under study. For unmixed ejecta, the temperature 
variation with time is that of Table 1 for each mass zone of the helium core whereas the initial 
gas number density at 100 days shows small changes due to the different initial molecular weight 
characterizing the chemical composition of each mass zone in the ejecta. 



4. RESULTS 

We now report on cluster abundances and masses for the 170 Mq and the 20 Mq fully 
microscopically-mixed ejecta in §4.1. We also present a critical assessment of the various exist- 
ing studies for fully-mixed ejecta to highlight some erroneous hypothesis assumed in certain models 
and the major problems encountered when using a CNT approach. Results for unmixed ejecta are 
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presented in §4.2. 



4.1. Fully-Mixed Ejecta 

As mentioned before, fully microscopically-mixed ejecta do not represent realistic cases of post- 



explosion SN environments as both theoretical explosion model s (IMiiller et al 



2003 



Hammer et al. 20091 ) and observations of SN remnants ( Douvion et al 



1991 



2001 



Kifonidis et al 



show evidence 

for the formation of inhomogeneities in the form of clumps and filaments in the ejecta. Because most 
of the existing studies on dust formation in SNe consider fully microscopically-mixed ejecta, and in 
order to compare and test the CNT and the kinetic approach, we present the results for such ejecta 
in Figure 2 for the 170 M© progenitor. The abundance variation with time for silica precursors 
and metal-bearing clusters are shown along with the abundances for the main atoms, ions, and 
secondary electrons. As already stressed in §2.1.1, silicon monoxide clusters (SiO) n will decompose 
at T > 900 K to |(Si02) n and g(Si) n clusters and can thus be considered as silica precursors. 
Until ~ 250 days, some atomic silicon is converted to SiO and SiS, building up a reservoir of 
SiO (see CD09 for more detail). For t > 250 days, the simultaneous and quick conversion of SiO 
to (SiO)2 and its further implication in the growth of larger clusters through the polymerization 
reaction [7] lead to the quick depletion of SiO and Si. The shape of the SiO abundance distribution 
is somewhat different than that presented in Figure 3 of CD09, and reflects the fact that in the 
present paper, the synthesis of silica precursors up to (SiO)s is followed whereas the growth up 
to (SiO)3 was studied in CD09. Silica precursors of smaller size than (SiO)s, including SiO, are 
converted to larger clusters but also replenished by the thermal fragmentation of (SiO)s reaction 
[8j thus keeping a bell-shaped abundance distribution until ~ 700 days when Si is depleted. Hence, 
via SiO formation, atomic Si is very efficiently and totally converted to (SiO)s at t > 700 days. 

Metal clusters form at late times, as illustrated in Figure 2b. Atomic iron, magnesium and 
aluminum are first destroyed by ionization to form Fe + , Mg + , and Al + ions before 200 days but 
recombine afterwards, and their abundance gradually reaches steady state, as opposed to atomic 
Si. The formation of pure magnesium and iron clusters is delayed to t > 700 days and is driven by 
the formation and destruction of (Mg)2 and (Fe)2 clusters, respectively (see reactions [271 Ell and 
[29]) . Indeed, the thermal deco mposition of (Mg)o and (Fe)2 is very efficient for temperatures in 



excess of 1000 K as derived bv lGiesen et al.l (|2003l ). but the decomposition reaction rate decreases 



dramatically for T < 1000 K, or t >700 days, leading to the build-up of pure magnesium and iron 
clusters as observed in Figure 2b. 

We also see that the 170 M@, fully-mixed ejecta does not synthesize metal oxide (FeO)4 and 
(MgO)4 clusters. The dominant formation process of FeO and MgO molecules is via reactions 
[9] and [I~3"l as O2 abundance sharply increases at t > 800 days. However, iron and magnesium 
oxidation competes with the fast (Fe)2,3,4 formation process reaction [27] and the formation of iron 
and magnesium oxides is therefore hampered. A similar competition between processes applies to 
the formation of metal sulphide clusters (FeS)4 and (MgS)4 clusters. The outcomes are very low 
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abundances for metal oxides and sulphides. As to alumina precursor A10, the destruction by He + 
precludes AlO to form at early times. However, He + abundances start decreasing at t ~ 500 days 
owing to the decrease with time of the ionization rate of He by Compton electrons (CD09). These 
lower He + abundances trigger the formation of O2 and the subsequent formation of AlO through 
reaction [T7] at t > 800 days. 

The derived cluster masses ejected at day 1000 are summarized in Table 4. (SiO)s clusters 
are by far the most abundant with a total mass of ~ 25.4 Mq, and decompose in 17.3 Mq of 
(SiO)5 and 8.1 Mq of (Si)s clusters, followed by small amounts of pure magnesium, (Mg)4, and 
iron, (Fe)4, clusters and aluminum oxide molecules, AlO. The total mass yield for clusters is ~ 25 
Mq, corresponding to 15% of the progenitor mass. We notice that no carbon chains and rings, 
and no silicon carbide (SiC)2 clusters form in this fully-mixed ejecta as atomic carbon is rapidly 
integrated into the synthesis of CO, and CS to a less extent (see CD09). 

Clusters abundances are shown in Figure 3 for the 20 Mq progenitor. Chemical processes 
similar to those at play in the ejecta of the 170 Mq progenitor are responsible for the nucleation 
of clusters. However, the initial chemical abundance of the ejecta at 100 days post-explosion is 
helium and carbon-rich, but oxygen and silicon-poor compared to that of the 170 Mq progenitor 
(see Table 2 from CD09 for detail). The high helium content impedes the early nucleation of silicon 
monoxide and pure metal clusters due to the large He + abundances in the ejecta. The formation 
of silicon monoxide and pure Si, Mg, and Fe clusters is thus postponed to t > 750 days when 
the He + abundance decreases. The lower gas density at these late times combined to lower initial 
mass fractions for Si and O result in smaller cluster masses synthesized in the ejecta, as illustrated 
in Table 4. Indeed, the total cluster mass equals ~ 0.16 M Q and corresponds to 0.8 % of the 
progenitor mass only. Thus, fully-mixed ejecta of primordial CCSNe are less efficient at nucleating 
solid clusters than PISNe due to their initial He-rich/metal-poor chemical composition. 

Prescriptions for the dust composition for the two mass progenitors under study are presented 
and discussed in 84.2.1. 



4.2. Testing Existing Studies 



The formation of dust at high redshift in S N fully-mixed ejec t a has been first modeled for 
zero-metallicity CCSNe by the pioneering study of lTodini &: Ferraral (|200ll . hereafter TF01). They 
described the condensation of dust grains using CNT and found that a 20 progenitor produced 



0.09 Mq of amorphous carbon dust. This formalism was further us ed by 



Schneider et al.1 



hereafter SFS04) to describe the formation of solids in PISNe, and by iBianchi &: Schneider! (120071 ) 
to study the p roduction of dust in local CCSNe and its survival to the passage of the reverse 
shock. Finally, IValiante et al.l (|2009l ) used TF01 dust results in a comparative study of dust main 
providers in the far and local universe. Another study on primordial SNe was carried out by NK03 
and considered both fully-mixed and unmixed ejecta of CCSNe and PISNe. 
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Firstly, as the ejecta composition, mixing, and physical parameters of NK03 are used in CD09 
and the present paper, we compared our results with those of NK03. Secondly, we test the validity 
of TF01 results in the framework of the chemical kinetic approach as TF01 results are used as a 
benchmark in many subsequent studies. Thirdly, we test the study of SFS04 in a similar context. 



4-2.1. The fully -mixed models by Nozawa et al. (2003) 

A direct comparison between results for our fully-mixed ejecta presented in §4.1 and those 
of NK03 is possible as NK03 ejecta models were used in CD09 and in the present study. NK03 
results for the fully-mixed ejecta of the 170 Mq and 20 Mq progenitors are listed in Table 5 
and compared to the present results. Despite identical initial conditions for the ejecta chemical 
composition and similar ejecta gas parameters, the present upper limits for the dust yields are 
lower than values derived by NK03, with a greater discrepancy (a factor of 2 to 5) for the 20 Mq 
progenitor. We also see that the chemical composition of the solids formed is totally different. The 
present models predict the formation of silica and pure silicon dust as the prominant condensates, 
along with alumina and pure magnesium and iron grains to a less extent. When full depletion of 
Fe and Mg is considered, fayalite and forsterite do also form. On the other hand, NK03 predict 
the formation of silica, forsterite/entastite (MgSiOa), magnetite (Fe304) and alumina. In their 
model, iron is trapped in magnetite whereas it is in the form of pure iron clusters in the present 
model. Also, NK03 do not form pure magnesium grains or pure silicon grains resulting from 
the disproportionation of silicon monoxide clusters at high temperatures. The only agreement we 
find regarding the dust composition is related to the lack of carbon and silicon carbide dust in the 
ejecta. NK03 assume that the whole atomic C is quickly transformed into CO, and we confirm from 
chemical kinetics that this is indeed the case for the fully microscopically mixed ejecta, whatever 
the initial mass of the progenitor. 

More differences arise in the condensation sequence derived by both studies. NK03 find for 
both cases that the condensation of alumina first proceeds, followed by forsterite/entastite, silica, 
and magnetite. For their 170 Mq progenitor, AI2O3 condenses at ~ 500 days, followed by Mg2Si04 
and MgSiOs at ~ 520 days, silica at ~ 550 days, and finally Fe304 at ~ 550 days. These various 
post-explosion times correspond to epochs when the nucleation rates of the various solids are 
maximum. However, as mentioned before, these nucleation rates are expressed as a function of 
physical quantities inappropriate to the description of molecular-size clusters, and the relevance of 
those condensation sequences have to be seriously questioned. Inspection of Figures 2 and 3 show 
different formation times for our cluster precursors, driven by chemical kinetics. In particular, the 
high initial silicon yields characterizing the 170 Mq progenitor foster a rapid formation of silicon 
monoxide clusters (hence silica and pure Si precursors) at high temperatures as soon a s 400 days, 



comparable to what is observed in ceramic- forming flames (jZachariah Sz Burguesslll994l ). The pure 
metallic clusters (Mg)4 and (Fe)4 on the other hand are forming at t > 700 days for the reasons 
explained in §4.1. Finally, AlO forms at day 900. Those strong discrepancies in formation times 
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reflect the fact that condensation from the gas phase is controlled by non-equilibrium chemically- 
controlled nucleation processes, each of which has its own temperature dependence. 

For their 20 M Q progenitor, NK03 find that Al 2 03 condenses first at ~ 380 days, followed by 
Mg 2 Si0 4 and MgSi0 3 at ~ 400 days, silica at ~ 420 days, and finally Fe 3 4 at ~ 430 days. This 
condensation sequence is similar to that for their PISN case though solid formation proceeds at 
slightly earlier times owing to the smaller ejecta temperatures. In the present study, nucleation 
times are totally different for the PISN and the CCSN cases, reflecting different chemistries at play 
in the ejecta. As already stressed in §4.1, the large He + content of the 20 Mq ejecta totally hampers 
cluster formation before day 700. Once again, these strong deviations are due to the inappropriate 
use of CNT applied to kinetically-commanded environments. 



4.2.2. The CCSN model by Todini & Ferrara (2001) 

We run our chemical network including the formation of molecules and nucleation of clusters 
for the initial elemental yields of the 20 Mq progenitor given by Table 2 of TF01 and their specific 
profiles for the temperature and number densit y of this ejecta. The eleme ntal composition used by 



TF01 is that derived from explosion models bv lWooslev &: Weaverl (j 19951 ) for their zero metallicity 



(Z = 0) progenitors. From TF01 Table 2, we notice for their 20 Mq (and 18 Mq) progenitor: 
(1) anomalously very low yields for heavy elements like Si, Mg, S and Al; (2) a carbon-to-oxygen 
ratio (C/O) larger than 1 compared to values for slightly lower or larger progenitor masses. Such 
very low mass yields for refractory elements and a C/O > 1 are not present in ot her explosion 



models like those of UN02 or the more recent calculations of lHeger Woosleyi (|2008l ) for a similar 
progenitor mass. Moreover, in those last two studies and for the specific 20 Mq progenitor mass, 
the elemental yields are found quite similar. Refractory elements like Si, Mg, or Al being all key 
players in the dust condensation processes, their paucity drastically changes the chemical nature 
of the dust that condenses. We also notice that the ejecta temperature profile used by TF01 is 
derived from comparison with the ejecta of SN1987A. The photospheric temperature is chosen equal 
to 5400 K and slowly decreases with time to reach 960 K at day 1000, as illustrated in Table 2. The 
ejecta temperature from 100 to 1000 days post-explosion derived by NK03 and used in the present 
study and in CD09 takes into account the deposition of energy by radioactive elements through the 
degradation of 7- and X-rays using a radiative transfer formalism coupled to the energy equation. 
As seen from Table 2, the resulting temperature variation with post-explosion time has a much 
sharper decrease with a higher initial temperature at day 100. As to gas number density, TF01 
mentioned a value of n = 1 X 10 8 cm -3 for temperatures less than 6000 K but do not provide a 
specific number density profile for their 20 Mq progenitor. We thus assume that the gas number 
density has the constant value of 10 s cm -3 over the timespan involved in our test study. 

The initial He/C-rich, Si/Mg-poor ejecta elemental composition and temperature profile used 
by TF01 imply that only carbon dust is expected to form if it can form. Indeed, they find that most 
of the dust formed for this specific 20 Mq progenitor is amorphous carbon with a AC dust mass yield 
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of ~ 0.09 Mq. When using a similar gas elemental composition, temperature and density profiles 
for our 20 Mq fully-mixed ejecta, our results show that the only molecules which are synthesized 
in the ejecta are C2 and CO, as illustrated in Figure 4. Their abundances with respect to the total 
gas number density are very low and comprised between 1 x 10~ 13 and 1 x 10 -11 from day 100 to 
day 1000 due to the overwhelming presence of He + . Indeed, the relatively high ejecta temperatures 
at day 300 and later fosters the ionization of He and impede He + recombination, hence sustaining 
a strong destruction pathway to molecular synthesis over the whole timespan under study. With 
the specific high helium yield in TF01 initial elemental composition, it is thus not surprising that 
molecule formation is hampered in this ejecta. As a consequence, no carbon chains, starting with 
the synthesis of C2, are formed and released at day 1000 owing to their destruction by He + , as 
seen in Figure 4. Our results are in total contradiction with TF01 carbon dust mass yield of ~ 
0.09 Mq. This deviation is due to the fact that their study considers the direct conversion of the 
initial carbon content decreased by its depletion in forming CO at steady state, without considering 
the ejecta full chemistry and the overwhelming destruction of small carbon clusters by He + . CD09 
show that a non-steady state chemistry commands the evolution of the ejecta gas phase until 1000 
days after explosion, and the present results illustrate again the importance of considering the 
ejecta chemistry and all chemical processes responsible for dust nucleation in any treatment of dust 
formation in SNe. 

In summary, we suggest that the high carbon dust yields derived by TF01 for fully-mixed 
ejecta of zero-metallicity 20 Mq CCSNe should be used with caution for they have been derived 
for questionable nucleosynthesis initial conditions and using the CNT approach. 



4.2.3. The PISN model by Schneider et al. (2004) 



Dust synthesis in Pop. Ill massive SNe is studied by SFS04. The CNT approach is used 
to describe the condensation of solids and SiC grains are not included among the condensates 
considered. The initial chemical compositions resulting from e xplosion nucleosyn t hesis for stellar 
progenitors with mass ranging from 140 to 260 Mq are those of lHeger &: Wooslevl (|2002l . hereafter 
HW02). Firstly, we notice that very low temperature profiles are obtained for the various massive 
progenitors considered, e.g., the 195 Mq progenitor ejecta has temperatures decreasing from ~ 
5400 K to ~ 50 K for 100 < t < 1000 days, a somewhat similar T variation than that found in 
CCSN ejecta. Secondly, inspection of SFS04 Figure 3 shows that the ejecta temperatures decrease 
with increasing progenitor masses. Low temperature profiles are questionable for PISNe, which 
are characterized by much larger explosion energies and 56 Ni mass synthesized in the explosion 
nucleosynthesis than CCSNe (HW02, UN02). As a matter of fact, solving for the radiative transfer 
equation taking into account the energy deposition by radioactive elements, NK03 find higher 
ejecta temperatures for PISNe than for CCSNe. Secondly, more massive PISNe are characterized 
by larger cent ral temperatures, explosion energies, and synthesized 56 Ni mass than their less massive 

, HW02, UN02), normally resulting in higher ejecta temperatures 
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due to the conversion of part of the explosion energy in thermal energy and the deposition of 
energy through 7-rays production from radioactive decay. As the condensation of solids in TF01 
and SFS04 formalism is assumed to start when the nucleation rates reach their maximum values, 
the low ejecta temperatures permit the condensation of grains as early as 150 days post-explosion 
where the gas densities are high enough to provide large final dust yields. 

Results from SFS04 and the present modeling are summarized in Table 6. SFS04 find that 
the amount of CO formed in the ejecta under steady state is equal to ~ 6 M and 0.01 M Q for 
the 170 and 260 M PISNe, respectively. Their dust mixture is dominated by silica, forsterite and 
amorphous carbon. They derive a total dust yield of ~ 39 M and ~ 70 M Q for the 170 M & and 260 
Mq progenitors, respectively. According to Table 6, the dust synthesized by a 260 M PISN also 
includes 8 Mq of magnetite Fe304, due to the large initial iron mass yield in the model of HW02. 
We apply our chemical kinetic approach to the SFS04 170 M and 260 M PISN ejecta models 
with similar initial compositions. Firstly, we find that a large molecular phase is synthesized in 
the ejecta for both models even when O2 depletion is considered in the derivation of the dust yield 
upper limits. For the 170 M progenitor, the molecular phase ranges from ~ 34 M and 51 M 
depending on the metal depletion, while it ranges from 33 Mq to 36 M for the 260 Mq CCSN. For 
both progenitors, molecules include CO, O2, and SO and very massive PISNe are less efficient at 
forming chemical species in their ejecta than less massive PISNe, owing to their hotter gas and their 
more important radioactivity at play, as already shown by CD09. The present dust mass upper 
limits are different from the dust mass yields derived by SFS04 for the two massive progenitors, 
although somehow similar in term of order of magnitude. However, strong disparities exist in both 
the chemical composition of the dust and the epoch of its formation. For example, SFS04 form 
amorphous carbon as the atomic carbon not trapped in CO turns into AC grains. They argue 
that some CO is destroyed by the 7-rays produced by 56 Co radioactive decay and produces this 
pool of free atomic carbon, with no quantitative assessment of this process in their ejecta. CD09 
show that this CO destruction pathway is of minor importance compared to alternative routes like 
neutral- neutral processes or destruction by He + . When kinetics is considered, the entire carbon 
initial content is quickly included into CO, resulting in much larger CO masses than those derived 
from steady-state, and the formation of small carbon chains is suppressed. No carbon dust can 
thus form in the present fully-mixed ejecta of PISN models as already seen in §4.1. 

As to other solids, we find that all the (SiO)s rings may be converted into silica, fayalite, 
forsterite, and pure silicon, or silica and pure silicon, dependent of the level of depletion of atomic 
Mg and Fe. Two other precursors are present in various amounts dependent of progenitor mass, AlO 
and (Fe)4, leading to the formation of alumina and pure iron grains. The large initial Fe elemental 
yield for the 260 Mq progenitor leads to very large amounts of pure iron dust at t = 1000 days when 
Fe depletion is not considered, as seen in Table 6. SFS04 find that silica, forsterite, amorphous 
carbon and magnetite are the dominant dust products for both progenitor masses. Their dust 
chemical composition is then very different from that derived in this study. Also important is the 
discrepancy in the condensation sequence. SFS04 find a condensation time window of 150 to 250 
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days, dependent of their progenitor masses. Nucleation epochs in the present study are controlled 
by chemical kinetics, reproducing dust synthesis in the laboratory. For example, the very low ejecta 
temperatures considered by SFS04 imply that the (SiO)s clusters form at 110 days post-explosion 
(or at temperatures close to those found in flames) whereas (Fe)4 are synthesized at t = 250 days 
and AlO forms at much later times (500 days) and lower temperatures. 

In summary, the SFS04 study uses very low ejecta temperatures for PISNe and the CNT 
approach to describe dust condensation. While the amount of dust synthesized is comparable to 
the upper limits we derive for similar progenitor masses, the present study points to a totally 
different composition for the condensed grains, in particular, the formation of pure iron grains and 
the lack of amorphous carbon and iron oxides like magnetite. Once again, caution should be taken 
when using CNT model outputs into star formation and chemical evolution studies in the early 
universe. 

4.3. Unmixed Ejecta 

We present in this section the results for the PISN 170 M Q and the CCSN 20 Mq unmixed 
ejecta. The stratification of the ejecta results from: (1) the nucleosynthesis of the explosion itself 
for the innermost mass layers; (2) successive nucleosynthetic phases occuring during the evolution 
of the stellar progenitor. Such a stratification pattern implies that each zone has its proper chemical 
composition and therefore its specific chemistry and chemical products. 

4.3.1. 170 M & progenitor 

Cluster abundances versus time for the 170 M & progenitor are presented in Figures 5 and 6 for 
Zones 1 and 2, respectively, while the dust mass yields ejected at 1000 days assuming Mg depletion 
or no depletion, as explained in §2.4, are summarized in Table 7. While Zone 1 is characterized by 
a unique Si/Fe-rich chemistry, Zone 2 contributes the most to the cluster budget from all zones. We 
see from Figure 5 that the first clusters to form at ~ 400 days are iron sulphide clusters, followed 
by pure silicon clusters at ~ 700 days and pure iron clusters at ~ 800 days. The FeS molecule 
is quickly formed from reaction [24] when its destruction occurs via thermal fragmentation and the 
reverse process of reaction [231 FeS formation leads to the depletion of molecular S2 at t > 300 days, 
as seen in CD09. The polymerization of FeS clusters occurs via reaction [26] at t > 400 days as 
thermal fragmentation remains important at early times due to the gas high temperatures. Once 
this process is switched off, the growth of the small cubic cluster (FeS)4 proceeds, with no other 
chemical destruction channels taking place. For pure silicon and iron clusters, formation of di-silicon 
and di-iron chains are assumed to proceed with a similar formation mechanism characterized by 
the three-body reaction [27] The formation of (Si)4 occurs at slightly earlier times than that of 
(Fe)4 owing to the large amounts of available atomic silicon not locked in SiS. The high Si initial 
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composition of Zone 1 spawns the large abundances of pure Si clusters, followed by iron sulphide 
and pure iron clusters. As to ions, the dominant species in Zone 1 is Fe + which forms from the 
charge exchange reaction between Si + and Fe, Si + being formed from the 7-rays ionization of atomic 
Si. Fe + is further destroyed by its recombination to Fe, which controls the electron population in 
the zone. 

In Zone 2, characterized by a high Si/O initial content, the conversion of atomic Si into SiO 
occurs for t < 250 days as seen in Figure 6, followed by the quick polymerization of SiO ring clusters 
at t > 350 days (equivalent to a gas temperature of ~ 3700 K). The further growth of (SiO)s keeps 
depleting the SiO content of the ejecta zone, resulting in large (SiO)s abundances. We stress here 
that the SiO distribution obtained in the present study and shown in Figure 6 is slightly different 
from that presented in CD09 for the same mass zone. This is due to the fact that we treat SiO 
ring polymerization up to 5 (SiO) units in this study whereas (SiO)3 was the end-product of cluster 
formation in CD09. The discrepancy results from the replenishment in SiO molecules during the 
destruction of (SiO^^s rings, leading to a higher SiO abundance distribution with time, although 
the SiO abundance trend remains the same. As to other condensates, AlO forms at T > 550 days 
from reaction 1171 once O2 has stopped to be efficiently destroyed by thermal fractionation. Finally, 
pure Mg clusters are synthesized at t > 800 days when their thermal fractionation has ceased. 

In table 7, we list the results on the upper limits for dust masses formed at 1000 days post- 
explosion for Zones 1-5 and for the two cases of total Mg depletion and no depletion. In this 
regard, the latter case gives almost directly the mass yields of the clusters synthesized at 1000 
days. As already stated, the dominant dust formation zones are Zones 2 and 3, followed by Zone 
1. From Table 7, we see that there exist two mass classes of condensates. The first class gathers 
the most abundant solids, all Si-bearing compounds, silica, forsterite, and pure silicon. Indeed, 
the prevalent condensate in both depletion cases is silica (Si02), forming essentially in Zones 2 
and 3. If Mg depletion is considered, forsterite (Mg2Si04) also forms in these zones. In both 
depletion cases, pure silicon clusters are too abundant owing to the disproportionation of silicon 
monoxide clusters, as explained in §2.1.1. This first class of solids represents ~ 99.5 % of the dust 
content. The second class of condensates gathers the minor compounds, alumina, iron sulphide, 
and pure magnesium and iron grains. All of those form separately from the first class and their 
mass yields are independent of the depletion case. This class accounts for ~ 0.5 % of the total dust 
budget. It is noteworthy that neither carbon grains nor periclase (MgO) and SiC precursors form 
in significant amounts in the unmixed ejecta of the 170 M & progenitor. Despite the presence of 
carbon in Zones 4 and 5, these zones are characterized by a C/O ratio less than unity, a very low 
atomic Si content, and by large amounts of helium in Zone 5. Although the formation of carbon 
chains is triggered by the rapid thermal conversion of CO into C2 at high temperatures and not 
by slow radiative association processes, as seen in CD09, those chains never form in significant 
quantities due to either their overwhelming oxidation by atomic O which hinders their growth in 
Zone 4, or the presence of He + in Zone 5. As to SiC precursors, they do not form owing to the 
small Si initial content of Zone 5 and the fact that the synthesis of gas-phase SiC competes with 
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the faster formation channels for C-chains, and suffers severe destruction by He + . 

Upper limits of the mass yields for each solid versus depletion case are shown in Figure 7 for the 
present 170 M & progenitor along with results from NK03 for a similar unmixed ejecta. According 
to Table 7, the total amount of dust formed for the 170 Mq progenitor ranges between 5.7 and 7.1 
Mq. This is well below the ~ 25 Mq of dust synthesized in NK03 170 Mq ejecta. As seen from 
Figure 7, our dust chemical composition also disagrees with that of NK03 who find that ~ 75 % 
of the grains are made of pure silicon, forsterite and silica grains, when the last 25 % are pure iron 
and iron sulphide grains, periclase, and amorphous carbon, the latter accounting for ~ 1% of the 
total dust mass. 



4-3.2. 20 M & progenitor 

Cluster abundances for Zones 1 and 2 of the 20 M Q progenitor are shown in Figures 8 and 9, 
respectively while upper limits for dust masses formed at 1000 days post-explosion for Zones 1-4 
and for the two cases of total Mg depletion and no depletion are presented in Table 8. The chemical 
processes at play in the formation of clusters are similar to those already highlighted in §4.3.1 for 
Zones 1 and 2 of the 170 Mq ejecta. In Figure 8 and for Zone 1, the prevalent clusters are, in 
order of decreasing abundances, (FeS)4, (Si)4, (SiO)3, and (Fe)4. A major departure arises from 
the large initial mass of oxygen compared to Zone 1 of the 170 Mq ejecta. In this Si-rich zone, the 
large O mass combined to the low gas temperature fosters the early nucleation of (SiO)3 clusters 
at t ~ 250 days, while SiO clusters are absent from Zone 1 of the 170 Mq ejecta. Moreover, Zone 1 
contains eight times more Fe than Zone 1 of the 170 Mq ejecta. All these discrepancies spawn the 
larger efficiency at making dust in Zone 1 for the CCSN, as shown in Table 8, where dust mass 
yields at 1000 days are presented. However, the formation of Fe-rich silicate like fayalite (Fe2SiC>4) 
is hampered in Zone 1 due to the fact that all the initial oxygen is trapped into the formation of 
SiO, leaving no oxygen for further oxidation by O2 or atomic O. Thus, the large amounts of (SiO)3 
result in silica formation. When Zone 2 forms essentially silica and/or forsterite, we see from Table 
8 that Zone 3 and Zone 4 do not contribute to the cluster budget. Zone 3 is Si-poor with a C/O 
ratio less than unity, implying that CO is the prevalent species to form in this zone (CD09). Zone 3 
is thus not conducive to cluster synthesis. Zone 4 is characterized by a large C/O ratio (29.5) and a 
high He content. These conditions preclude the synthesis of C-chains and rings, and the subsequent 
AC grain condensation because of the destruction of small clusters by He + , as explained in §4.4.2. 
For SiC precursors, conclusions similar to those drawn for the 170 Mq case apply, and no silicon 
carbide dust forms in Zones 3 and 4. 

Upper limits of the mass yields for each solid versus depletion case are shown in Figure 10 
along with the results from NK03 for a similar unmixed ejecta. According to Table 8, the total 
amount of dust formed for the 20 M Q progenitor ranges between 0.103 and 0.154 Mq. This is well 
below the ~ 0.58 Mq of dust synthesized in NK03 170 Mq ejecta. As seen from Figure 10, our dust 
chemical composition disagrees with that of with NK03 who find that ~ 64 % of the dust mass is 
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in the form of pure silicon, forsterite and silica grains, 16 % and 5 % are iron sulphide and pure 
iron grains, respectively, and 10 % is in the form of amorphous carbon. 



4-3.3. Depletion efficiencies 

Elements initially in atomic form are first depleted in molecules and molecular clusters as a 
result of nucleation chemistry. To assess this trapping independent of the depletion cases defined 
in §2.5, we define the efficiency x as 

X = Nciuster(y)/N To tal(y), (41) 

where y is the element, Nci uster (y) is the total number of y atoms locked up in one type of cluster, 
and NTotai{y) is the initial total number of y atoms in the ejecta. The efficiency x measures how 
the chemistry locks up elements in the nucleation processes of molecular clusters. It is not linked to 
the full depletion of Mg and Fe during the condensation phase of dust seeds. Depletion efficiencies 
are shown as a function of zoning in Table 9 for the two progenitor masses. 

For both progenitor masses, the depletion of Si atoms into silicon oxide precursors leading to 
silica and pure silicon clusters is efficient in silicate forming zones and starts at early times when 
Si is converted into SiO and small SiO clusters (see Figures 6 and 9). The depletion is then total 
at day 1000 as seen in Table 9 and §4.3.1 and §4.3.2. This is in agreement with the steeper fading 
of the 1.65 /iin [S il] line obse r ved in SN1987A compared to the 4571 A [Mgl] and the 6300 A [OI] 
lines observed by Lucy et al. (jl989l ). In the Fe-rich, O-poor Zone 1, the conversion of Si into pure 



silicon clusters proceeds with a low efficiency of ~ 1% to 3% due to the prevalent conversion of Si 
in SiS molecules with an efficiency ranging from 26 to 28%. 

Aluminum also suffers total depletion into molecular AlO. The formation of alumina clusters 
has not been modele d here but AlO should react with itself and O2 to build up small AI2O3 clusters 



(jCatoire et al.ll2003l ). A total depletion of Al into alumina clusters is thus to be expected. 



Carbon is massively depleted into CO in the zones where it is present. In the C-rich Zones 4 
and 5 of the 20 and 170 Mq ejecta, respectively, its depletion depends on the presence of He + , as 
fully explained in §4.4.2. In Zone 5 of the 170 Mq PISN, which has a C/O ratio of 0.66, most of 
the carbon is tied up in CO once the He + abundance decreases for t > 600 days (see Figure 7 of 
CD09). Conversely, in Zone 4 of the 20 Mq CCSN characterized by a C/O ratio of 29.5, C remains 
in atomic form due to the overabundance of He, hence He + . However, atomic carbon is almost 
entirely depleted in the carbon ring C10 if He is not microscopically mixed to the gas, as seen in 
§4.4.2. 

Atomic oxygen is depleted in molecular O2 with efficiency values ranging from 12% to 61% 
in all zones where O is present, and into CO (x = 29 — 43%) where atomic carbon is abundant. 
The formation of silica clusters locks up atomic oxygen to efficiency values spanning ~ 3% to 10% 
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only. Thus part of the decline in the oxygen lines is due to molecular formation rather than dust 
condensation. 

Iron is depleted in the formation of FeS clusters in Zone 1 of both ejecta. Its depletion 
efficiency in FeS reaches the high level of 44% for the 20 Mq progenitor ejecta. We also see that Fe 
is never efficiently transformed into Fe clusters. To assess if this low efficiency in nucleating iron 
small clusters is due to the presence of sulphur and the formation of FeS, we ran the chemistry 
for the S-free Zone 1 of the 20 Mq progenitor. We find that the depletion of iron into small Fe 
clusters stays very low (~ 0.3%) and that most of the iron remains in atomic form. As to Mg, 
it is seldom depleted by the nucleation of small clusters except in Zones 3 and 2 for the 170 and 
20 M & ejecta, respectively, with depletion efficiencies of 1 and 2%. However, the coagulation into 
large silica clusters may trapped atomic Mg or Fe into the cluster lattices to form forsterite and 
fayalite, respectively, a situation described by our 'total depletion' case defined in §2.1. However, 
this depletion processes are beyond the aims of the present paper and will be presented in a 
forthcoming publication. 

4.4. Critical Parameters under Study 

We now explore the impact of two critical parameters in the cluster nucleation process, i.e., 
the nucleation rate of reactions [7] for silica clusters and the impact of He diffusion on carbon dust 
formation in the outer zones of the ejecta. 

4-4-1- Pressure Dependence of Nucleation Reactions 

A critical parameter to the formation of silicates and metal oxides is the rate at which small 
clusters react one to another to nucleate small dust seeds that can later coagulate and condense 
into lager dust grains. In previous chapters, we assume that reactions [7J Q21 [161 G31 and [26] all 
proceed with similar rates listed in Table 2, as the only information available on nucleation rates 
is for small silicon oxide clusters (ZT93). To test the impact of varying the nucleation rates of 
small clusters of silicon oxide, metal oxides and sulphides, we consider a case study describing the 
nucleation of small silicon oxide clusters in Zone 2 of the 170 Mq progenitor unmixed ejecta. The 
rates for reactions such as reactions [7J are listed in Table 2 and are taken from the ZT93 study for 
the 10 -2 atm pressure case. Using this series of rates as benchmarks, we now study the impact of 
increasing and decreasing by a factor 100 the benchmark rates. The increase by 100 corresponds to 
the 1 atm pressure case of ZT93 and we calculate the abundances normalized to total gas number 
density of clusters as a function of post-explosion time. 

Abundances for the (SiO)s clusters and gas-phase SiO are presented in Figure 11 for the 
different rates of reactions [71 As the nucleation rate decreases (e.g., in the 10 -4 case for which all 
rates have been decreased by a factor 10~ 2 compared to our benchmark case), the formation of 
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clusters is postponed to late times while the SiO abundance increases first and starts dropping once 
(SiO) n clusters form. However, the abundance of the (SiO)s cluster is unchanged at day 1000. Not 
shown in Figure 11, the O2 and CO abundances follow similar trends, implying that the lower the 
rates, the larger the molecular phase and its mass in the ejecta while the final cluster mass yields 
at day 1000 are similar for all cases. As seen in §4.3.3, elemental Si and SiO are totally depleted in 
the process of silica and silicate dust formation. Therefore, observations of Si and SiO optical and 
IR lines in local young SNe are good indicators of the timing and the efficiency of the nucleation 
of silicate dust in SN ejecta, once excitation effects are taken into account in the analysis of line 
declines as a function of time. 



How Can Carbon Dust Form? 



No carbon rings form in the unmixed ejecta of the 170 and 20 Mq progenitors owing to the 
large amounts of He and hence He + in C-rich Zones 5 and 4, respectively. However, evidence 
for the presence of carbon grains in CCSNe ejecta exists from various sources of investigation . 
Graphite grains bearing the isotopic anomalies of CCSNe are found in meteorites (jZinnerll2006l ). 
Furthermore, models aiming at explaining the IR excess emission of several low-z CCSNe and SN 
remnants show that a morphous carbon is a necessary component of the dust mix ture required 
to reproduce IR data (jSugerman et al.l 12000 : lErcolano et al.1 120071 : iRho et al.l 120081 ). In order to 
assess the impact of He/He + on the formation of carbon chains and rings, we re-investigate the 
chemistry of Zone 4 of our 20 Mq progenitor, which is characterized by a C/O ratio of 29.47. 
Keeping the initial metal mass constant, we gradually increase the amount of He present in the 
zone up to its value assumed in the previous section (0.884 Mq according to Table 3 of CD09). 
We consider the following increase factors: 0% (= no He), 1%, 20%, 40%, and finally 100% of the 
He value of Zone 4. For increasing He masses, the initial total number of elements in the gas also 
increases while the partial number of metals is kept constant. Doing so, we aim at reproducing the 
chemical composition of a C-rich clump in the ejecta partially embedded in a helium mushroom 



cap, a configuration commonly found in 2 D/3D modeling of CCSN explosion (jMiiller et al. 



1991 



Kifonidis et al.l 120031 : lHammer et al.l |2009j) . We then allow some diffusion of He inside the C-rich 
Zone 4. 

The abundance of the ring C10 as a function of post-explosion time for our various He diffusion 
values are shown in Figure 12. We see that the larger the He content of the zone, the later the 
formation and the smaller the abundances of C10 rings. Indeed, for very large He content, the 
amount of He + in the zone is large (x(He + ) ~ 0.1% x 2; (He), where x(He) and x(He + ) are the 
abundances of He and He + , respectively). The destruction of small carbon chains like C2 then 
proceeds until He + abundances diminishes at late times owing to a decrease in He ionization by 
Compton electrons and the recombination of He + . The resulting chain and ring formation is thus 
delayed. When the zone is free of He, the rapid conversion of the atomic carbon not locked in 
CO into C2 via the radiative association reaction [30] triggers the formation of larger chains, with 
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instantaneous closure of the end-product ring Cio- For the He- free Zone 4, the final abundance of 
Cio with respect to total gas number density at 1000 days is 0.55, equivalent to a AC dust mass 
yield of 0.0145 Mq in which ~ 95% of the initial carbon is depleted. When He diffusion occurs, 
the conversion of free atomic carbon into Cio is as effective as for the He- free case, with similar 
final mass yield and formation efficiency for the Cio ring, except that the conversion takes place at 
later times until the He + abundance decreases. 

In summary, the build-up of carbon clusters in SN ejecta is significantly hampered by oxidation 
in zones where the C/O ratio is less than unity, and by the microscopic mixing of He + in the carbon- 
rich zones of the ejecta. More generally, He + ha s a detrimental impact on molecular formation 



in SN ejecta as first shown by iLepp et al.l (|l990l ) and CD09. In C-rich, He-free zones, carbon 



condenses with a high efficiency. Any formation of carbon dust in SNe must thus occur in C-rich 
inhomogeneities deprived of He + and is thus highly dependent of the mixing induced by the SN 
explosion. 



5. SUMMARY AND DISCUSSION 

We have investigated the synthesis of molecular clusters precursors to dust in Pop. Ill supernova 
ejecta and show that dust clusters form despite the harsh, post-explosion environments. Of special 
interest are the following points: 

• In general terms, we outline that a classical nucleation theory (CNT) formalism is not appro- 
priate to the descrition of dust formation in SNe ejecta as it involves hypothesis which do not 
hold for SN ejecta out of chemical equilibrium and concepts not prevailing for small molecular 
clusters. We highlight some specific problems encountered by previous stu dies of fully-mixed 



SN ej ecta using CNT, specifically: (1) erroneous initial conditions for the lTodini &: Ferrara 



(|200ll ) 18 and 20 Mq progenitor models leadi ng to the format i on of carbon dust in Pop. Ill 



SNe; (2) very low ejecta temperatures used bv lSchneider et al.l (|2004! ) for their PISN models 



A chemical kinetic approach applied to the 20 Mq TF01 model makes the synthesis of dust 
merely impossible because of the large He content of the ejecta. When applied to SFS04 
models, the chemical kinetic approach points to a totally different chemical composition for 
the dust than that derived by SFS04. In both cases, chemical kinetics prevents the formation 
of carbon dust. 

Fully microscopically- mixed and unmixed models produce smaller dust masses ejected at day 
1000 after explosion than those de r ived by studies using CNT and similar ejecta conditions 
such as the study by Nozawa et al. ( 20031 ). Dust mass upper limits ranging from 25 to 33 Mq 



and 0.16 to 0.33 Mq are derived for the fully-mixed ejecta of 170 Mq and 20 M progenitors, 
respectively. These upper limits depend on the level of Fe/Mg depletion and lie between 16% 
and 21 % of the 170 Mq progenitor mass and 0.8% and 2% of the 20 Mq progenitor mass. 
NK03 find that 170 and 20 Mq progenitors form dust mass yields representing 21% and 4% 
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of the progenitor mass, respectively. Our upper limits are thus down by a factor of 2 to 5 for 
the 20 Mq progenitor. Our dust composition is dominated by silica (and silicates when metal 
depletion is considered), pure silicon, alumina and pure iron. No iron oxides, amorphous 
carbon dust, or silicon carbides form. These dust compositions disagree with those derived 
by NK03. 

• For the unmixed ejecta, we find dust upper limits dependent of the degree of metal depletion 
ranging from 5.7 to 7 Mq for the 170 Mq progenitor, and from 0.09 to 0.14 Mq for the 20 Mq 
progenitors. These limits are much smaller than the dust masses found by NK03 by factors 
of ~ 4 and ~ 5 for the 170 and 20 Mq progenitors, respectively. The chemical composition 
of the dust is dominated by silica/silicates, pure silicon dust and iron sulphide grains. No 
magnesium and ferrous oxides, pure iron grains, silicon carbides, and amorphous carbon dust 
form, as opposed to the results from NK03. 

• Negligible amounts of carbon dust form in unmixed ejecta, despite the large carbon mass 
initially present in the outer ejecta zone of the 20 Mq progenitor. We show that the presence 
of He + hampers the formation via small carbon chains and rings even for high C/O ratios. 
Consequently, carbon dust formation in SNe must occur in inhomogeneities where He/He + is 
not microscopically mixed with the gas. In zones where the C/O ratio is less than unity, ther- 
mal fragmentation of CO triggers the formation of C2 and longer carbon chains, as opposed 
to radiative association reactions, but the oxidation of chains prevent those from growing 
into rings. No silicon carbide precursors form in mixed and unmixed ejecta as their synthesis 
channels compete with those for carbon precursors. Furthermore, the extremely low silicon 
content of the outer carbon-rich layers precludes SiC nucleation in unmixed ejecta. 

• The depletion of metals resulting from the formation of molecules and the nucleation of dust 
clusters prior to any depletion stemming from the dust condensation phase varies according 
to the elements and their local environment. In a O-rich gas, Si is depleted at 100% in the 
formation of silicon monoxide clusters, when a S-rich/O-poor gas depletes Si to a 30% level, 
mainly in gas phase SiS. Al is also totally depleted by the formation of gas-phase aluminum 
oxide. However, Fe is scarcely depleted except in the innermost mass zone of the CCSN ejecta 
where its conversion to FeS clusters accounts for an efficiency of 40%. Mg does not suffer 
severe nucleation depletion and gets depleted into small pure Mg clusters to a level of ~ 2 %. 

• The molecular content of the ejecta before 1000 days after explosion is highly dependent 
of the nucleation reaction rates at which small clusters grow. Smaller rates imply a larger 
molecular content for the ejecta. From a chemical point of view, optical and IR lines of [Sil] 
and SiO are thus direct tracers of the timing and efficiency of silicate condensation in local 
young SNe ejecta. 

In spite of the fact that the amounts of dust clusters forming in our unmixed PISN and CCSN 
ejecta are smaller than those derived by existing studies, the upper limits for dust are still high 



- 30 - 



compared to what is observed in local CCSNe. Observations with Spitzer (jSugerman et al.l 120061 ; 
Kotak et al.l 120091 ) point to dust mass yields not exceeding 10~ 2 Mq i.e., a factor 10 down to 
the upper limit we derived for a primordial 20 Mq CCSN. A possible reason for this discrepancy 
is the existence of cool dust escaping IR detection and its contribution to the final dust budget 
([Gomez et al.ll2009l ). Another possibility is that existing studies and the present work do not prop- 
erly consider real SN ejecta for their dust synthesis models. SN ejecta and remnants are known 
to be highly inhomogeneous, with the presence of clumps, knots and filaments forming after the 
explosion as a result of Rayleigh- Taylor instabilities developing at the interfaces of the nucleosyn- 
thesis zones of the progenitor. On the one hand, fully microscopically-mixed ejecta are highly 
unrealistic since the mixing going inward to the deep core layers as a result of Rayleigh- Taylor 
instabilities is macroscopic. On the other hand, unmixed ejecta do not describe the mixing with 
(non)-adjacent layers. Such a microscopic mixing is necessary to explain the formation of silicon 
carbide an d silicon nitride pre-solar inclusions bearing SN isotopic signatures found in meteorites 
(|Ottll2003l ). Therefore, we see that mixing in knots and clumps is a key parameter in the determina- 
tion of the chemical composition of the dust and in its formation efficiency. This is best exemplified 
with amorphous carbon, as we see that carbon chains, rings and grains can only grow in primor- 
dial CCSNe if the formation locus is carbon-rich and deprived of He + . More generally, any He + 
microscopic mixing impe des the format i on of molecules and molecular clusters, as seen in CD09. 
The recent 3D model bv lHammer et al.l (|2009l ) shows the formation of O-rich blobs encompassing 
up to 10 -2 Mq a few days after explosion. Depending on the blob initial chemical composition 
resulting from this mixing, a certain type of dust will form during the blob expansion. We see that 
any atomic Si will be totally depleted in a O-rich inhomogeneity due to the formation process of 
silica/silicate clusters, and the final amount of condensates will essentially depends on the initial 
Si/O ratio characterizing the inhomogeneity. If part of the Si resides in O-free blobs, its efficiency 
at forming pure Si clusters is rather low (~ 2% according to Table 7), and silicon will then stays 
in atomic form or may be depleted in SiS to a level of ~ 27 % if sulphur is present. Thus, mixing 
itself by the way of the chemical processes involved in the formation of molecules and dust clusters 
in inhomogeneities contributes towards the limitation of the total amount of dust synthesized in 
PISN and CCSN ejecta. 

Pop. Ill stars exploding as PISNe are definitely the first dust makers in the early universe 
and the first dust enrichment of the primordial gas should be silicon/silica/silicate-rich to a level 
of 99 % of the total dust mass ejected. No carbon grains are formed as atomic carbon is trapped 
in molecular CO in the outer layers of the PISN ejecta. Therefore, carbon dust at high redshift 
must come from different stellar sources than PISNe and its formation may be delayed until the 
lower-mass Po p.II.5/II stars have evolved to their AG B stage or as Supergiants, depending on their 



initial masses (Dwek et al. 



2007 



Valiante et al.ll2009i ) . PIS N dust will be further reprocessed by the 



reverse shock several thousands years after the explosion. iNozawa et al.l (|2007l ) find for an initial 
170 Mq progenitor a grain survival rate ranging from ~ 0.02% to 48% of the initial dust mass 
depending on the pre-shock gas number density and the initial explosion energy. Thus between 
~ 0.001 to 3 Mq of Si-based dust will be injected into the dense shell of the radiative phase of 
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a 170 M Q PISN explosion. This dust will further provide cooling for the gas fragmentation and 
Pop. II. 5 star formation. As an extension of the present work, dust mass yields, composition, and 
grain size distributions will be derived using a Monte-Carlo formalism for dust condensation in 
PISN ejecta. The impact on cooling of this Si-based dust and its related molecular phase as well 
as its contribution to Pop. II. 5 star formation will be further explored. 

We thank the anonymous referee for comments that have contributed to the improvement of 
the manuscript. IC acknowledges support from the Swiss National Science Foundation through the 
Maria Heim-Vogtlin subsidies PMPD2-114347 and PMPDP2-1241159. 
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Table 1. Molecular clusters and their related dust products considered in this study. 



Parent molecule 


Clusters 


Dust 


Symbol 


SiO 


(SiO) 2 to (SiO) 5 


Silica 


SiOa 






Forsterite 


Mg 2 Si0 4 a 






Fayalitc 


Fe 2 Si0 4 a 


MgO 


(MgO) 2 to (MgO) 4 


Magnesium oxide or Periclase 


MgO 


FeO 


(FeO) 2 to (FeO) 4 


Iron oxide or Troilitc 


FeO 


AlO 




Alumina 


AI2O3 


FeS 


(FeS) 2 to (FeS) 4 


Solid iron sulphide 


FcS 


MgS 


(MgS) 2 to (MgS) 4 


Solid magnesium sulphide 


MgS 


Fe 


(Fc) 2 to (Fe) 4 


Solid iron 


Fe 


Mg 


(Mg) 2 to (Mg) 4 


Solid magnesium 


Mg 


Si 


(Si) 2 to (Si) 4 


Solid silicon 


Si 


C 


C 2 to Cio 


Amorphous carbon 


AC 


SiC 


(SiC) 2 


Silicon carbide 


SiC 



a Formcd in the full metal (Mg and/or Fe) depletion case. 
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Table 2. Reaction rates for the chemical processes involved in the nucleation of molecular 
clusters a . Other reactions included in the total network are listed in Cherchneff & Dwek (2009). 



Chemical processes 


Aij 


V 




Reference b 


Aluminum - Al 



Al 

A2 

A3 

A4 

A5 

A6 

A7 

A8 

A9 

A10 

All 



Al + 2 
Al + C0 2 
Al + H 2 
Al + O + M 
AlO + O 
A10+ He+ 
AlO + He+ 
AlO + M 

AlO + e C ompton 

AlO + eZ 

Compton 

AlO + e^ 

Compton 



AlO + CO 
AlO + H 2 
AlO + M 
Al + 02 
Al + 0+ + He 
Al+O +He 
Al + O + M 

A' + O + ^Compton 

A1+ + O + c" 

Compton 

Al + 0+ + c" 

| Compton 



1.61X10- 10 -1 NIST 

2.89XKT 10 3221.0 Garland ct al. 1992 

1.59xKr 10 2868.6 Mc Lean ct al. 1993 

8.27xl0" 31 NIST 

1.93X10" 11 -0.5 1199.9 Garland ct al. 1992 

8.60xl0" 10 NIST 

8.60xl0" 10 NIST 

4.40xl0" 10 98600.0 Same as S2 

5.81xl0" 6 3464.1 See CD09 c 

2.95xl0" 6 3464.1 See CD09 c 

9.48xl0" 7 3464.1 See CD09 c 



Silicon - Si 



Si 2 2 


4.61x10" 


17 





-2821.4 


Zachariah & Tsang 1993 


SiO + SiO + M 


4.40x10" 


10 





98600.0 


Reaction TF15 in CD09 


S13O3 


2.24x10" 


15 





-2878.9 


Zachariah & Tsang 1993 


SiO + Si 2 2 + M 


4.40x10" 


10 





98600.0 


Zachariah & Tsang 1993 


S14O4 


1.53x10" 


14 





-2386.8 


Zachariah & Tsang 1993 


S14O4 


1.53x10" 


14 





-2386.8 


Zachariah & Tsang 1993 


S13O3 + SiO 


1.53x10" 


14 





-2386.8 


Zachariah & Tsang 1993 


S13O3 + SiO 


4.40x10" 


10 





98600.0 


Same as S2 


Si 2 2 + Si 2 O a 


4.40x10" 


10 





98600.0 


Same as S2 


Si 5 5 


1.53x10" 


14 





-2386.8 


Same as S5 


Si 5 5 


1.53x10" 


14 





-2386.8 


Same as S5 


Si 4 4 +SiO 


1.53x10" 


14 





-2386.8 


Same as S5 


Si 2 2 + Si 3 3 


4.40x10" 


10 





98600.0 


Same as S2 


SiO + Si 4 4 


4.40x10" 


10 





98600.0 


Same as S2 


Si 2 +M 


2.76x10" 


29 





-2821.4 


Same as F23 


Si + Si +M 


7.14x10 


-5 





17800 


Same as F24 


Si 3 +M 


2.76x10" 


29 





-2821.4 


Same as F25 


Si + Si 2 +M 


1.66x10 


-5 





19200 


Same as F26 


Si 4 +M 


2.76x10" 


29 





-2821.4 


Same as F27 


Si + Si 3 +M 


8.30x10 


-7 





21600 


Same as F28 


Si 3 +Si 


8.30x10" 


10 








Same as F29 


Si 4 


8.30x10" 


10 








Same as F30 


SiC + 


5.43x10 


-7 


-1.50 


57200 


Same as NN57 in CD09 


Si + CO 


5.99x10" 


10 


0.0 


1420 


Same as C38 


SiC + 


5.43x10 


-7 


-1.50 


57200 


Same as S23 


SiO + c 


5.99x10" 


10 


0.0 


1420 


Same as C38 


SiC + c 


5.99x10" 


10 


0.0 


1420 


Same as C38 


Si 2 C 2 


4.61x10" 


17 





-2821.4 


Same as SI 


Si 2 C 2 + M 


2.48x10" 


21 


-3.0 





Same as C37 


SiC + SiC+M 


4.40x10" 


10 





98600.0 


Same as S2 


Iron 


- Fe 










FeO+ O 


2.09x10" 


10 





10199.6 


Akhmadov ct al. 1988 


FeO + CO 


5.38x10" 


10 





15033.1 


NIST 



51 SiO + SiO 

52 Si 2 2 + M 

53 SiO + Si 2 2 + M 

54 Si 3 3 + M 

55 SiO + Si 3 3 + M 

56 Si 2 2 + Si 2 2 

57 Si 2 2 + Si 2 2 

58 Si 4 4 + M 

59 Si 4 4 + M 

510 SiO + Si 4 4 + M 

511 Si 2 2 + Si 3 3 + M 

512 Si 2 2 + Si 3 3 + M 

513 Si 5 5 + M 

514 Si 5 5 + M 

515 Si + Si + M 

516 Si 2 + M 

517 Si + Si 2 + M 

518 Si 3 + M 

519 Si + Si 3 + M 

520 Si 4 + M 

521 Si 2 + Si 2 + M 

522 Si 2 + Si 2 + M 

523 Si + CO 

524 0+ SiC 

525 C + SiO 

526 0+ SiC 

527 Si + C 2 

528 SiC + SiC 

529 SiC + SiC + M 

530 Si 2 C 2 + M 



Fl 
F2 



Fe+ 2 
Fc + C0 2 
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Table 2 — Continued 



Chemical processes 



Aij 



Reference 



F3 


FeO + O 


F4 


FeO + FeO 


F5 


Fe 2 2 + M 


F6 


FeO + Fc 2 2 


F7 


Fe 3 3 + M 


F8 


FeO + Fe 3 3 


F9 


Fe 4 4 + M 


F10 


Fe 4 4 + M 


Fll 


Fe+ SO 


F12 


Fe+ S 2 


F13 


FeS + O 


F14 


FeS + S 


F15 


FeS + FeS 


F16 


Fe 2 S 2 + M 


F17 


FeS + Fc 2 S 2 


F18 


Fe 3 S 3 + M 


F19 


FeS + Fe 3 S 3 + M 


F20 


Fe 4 S 4 + M 


F21 


Fe 4 S 4 + M 


F22 


Fe + Fe + M 


F23 


Fe 2 + M 


F24 


Fe + Fe 2 + M 


F25 


Fe 3 + M 


F26 


Fe + Fe 3 + M 


F27 


Fe 4 + M 


F28 


Fe 2 + Fe 2 + M 


F29 


Fe 2 + Fe 2 + M 



Ml 


Mg+ 2 


M2 


Mg + C0 2 


M3 


MgO + O 


M4 


MgO + MgO 


M5 


Mg 2 2 + M 


M6 


MgO + Mg 2 2 


M7 


Mg 3 3 + M 


M8 


MgO + Mg 3 3 + M 


M9 


Mg 4 4 + M 


M10 


Mg 4 4 + M 


Mil 


Mg+ SO 


M12 


Mg+ S 2 


M13 


MgS + O 


M14 


MgS + S 


M15 


MgS + MgS 


M16 


Mg 2 S 2 + M 


M17 


MgS + Mg 2 S 2 


M18 


Mg 3 S 3 + M 



Fe + 2 


4.60x10" 


10 





360 


Self & Plane 2003 


Fe 2 2 


4.61x10- 


17 





-2821.4 


Same as SI 


FeO + FeO + M 


4.40x10" 


10 





98600.0 


Same as S2 


Fe 3 3 


2.24x10" 


15 





-2878.9 


Same as S3 


FeO + Fc 2 2 + M 


4.40x10" 


10 





98600.0 


Same as S4 


Fe 4 4 


1.20x10- 


32 





2160.0 


Same as S5 


FeO + Fc 3 3 + M 


4.40x10- 


10 





98600.0 


Same as S8 


Fc 2 2 + Fe 2 2 + M 


4.40x10- 


10 





98600.0 


Same as S9 


FeS+ O 


2.09x10" 


10 





10199.6 


Same as Fl 


FeS+S 


2.09x10- 


10 





10199.6 


Same as Fl 


Fe + S0 2 


1.00x10- 


■10 








E 


Fe + S 2 


1.00x10- 


■ 10 








E 


Fc 2 S 2 


4.61x10- 


■17 





-2821.4 


Same as F5 


FeS + FeS + M 


4.40x10- 


10 





98600.0 


Same as F6 


Fe 3 S 3 


2.24x10- 


■15 





-2878.9 


Same as F7 


FeS + Fe 2 S 2 + M 


4.40x10- 


■10 





98600.0 


Same as F8 


Fe 4 S 4 


1.20x10- 


-32 





2160.0 


Same as F9 


FeS + Fe 3 S 3 + M 


4.40x10- 


■10 





98600.0 


Same as F10 


Fc 2 S 2 + Fc 2 S 2 + M 


4.40x10- 


■10 





98600.0 


Same as Fll 


Fc 2 +M 


2.76x10- 


■29 





-2821.4 


Giesen ct al. 2003 


Fe + Fc +M 


7.14x10 


— 5 





17800 


Giesen ct al. 2003 


Fc 3 +M 


2.76x10- 


29 





-2821.4 


Giesen ct al. 2003 


Fe + Fe 2 +M 


1.66x10 


—5 





19200 


Giesen ct al. 2003 


Fe 4 +M 


2.76x10- 


29 





-2821.4 


Giesen ct al. 2003 


Fe + Fe 3 +M 


8.30x10 


— 7 





21600 


Giesen et al. 2003 


Fe 3 +Fc 


8.30x10- 


■10 








Giesen et al. 2003 


Fe 4 


8.30x10- 


■10 








Giesen et al. 2003 


Magnesium 


- Mg 










MgO+ O 


1.18x10- 


10 





5769.1 


Kashireninov et al. 1982 


MgO + CO 


5.38x10- 


10 





15033.1 


Same as F2 


Mg + 2 


1.00x10- 


11 








E 


Mg 2 2 


4.61x10- 


17 





-2821.4 


Same as SI 


MgO + MgO + M 


4.40x10- 


10 





98600.0 


Same as S2 


Mg 3 3 


2.24x10- 


15 





-2878.9 


Same as S3 


MgO + Mg 2 2 + M 


4.40x10- 


10 





98600.0 


Same as M5 


Mg 4 4 


1.20x10- 


32 





2160.0 


Same as S5 


MgO + Mg 3 3 + M 


4.40x10- 


■10 





98600.0 


Same as S8 


Mg 2 2 + Mg 2 2 + M 


4.40x10- 


10 





98600.0 


Same as S9 


MgS+ O 


2.09x10- 


10 





10199.6 


Same as Fll 


MgS+S 


2.09x10- 


10 





10199.6 


Same as F12 


Mg + S0 2 


1.00x10- 


10 








E 


Mg + S 2 


1.00x10- 


10 








E 


Mg 2 S 2 


4.61x10- 


17 





-2821.4 


Same as F15 


MgS + MgS + M 


4.40x10- 


10 





98600.0 


Same as F16 


Mg 3 S 3 


2.24x10- 


15 





-2878.9 


Same as F17 


MgS + Mg 2 S 2 + M 


4.40x10- 


■10 





98600.0 


Same as F18 
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Table 2 — Continued 



Chemical processes 



Aij 



Reference 



M19 MgS + Mg 3 S 3 + M 

M20 Mg 4 S 4 + M 

M21 Mg 4 S 4 + M 

M22 Mg + Mg + M 

M23 Mg 2 + M 

M24 Mg + Mg 2 + M 

M25 Mg 3 + M 

M26 Mg + Mg 3 + M 

M27 Mg 4 + M 

M28 Mg 2 + Mg 2 + M 

M29 Mg 2 + Mg 2 + M 



Mg 4 S 4 

MgS + Mg 3 S 3 
Mg 2 S 2 + Mg 2 S 
Mg 2 +M 



M 
■2 + M 



Mg - 
Mg 3 
Mg - 
Mg 4 
Mg - 
Mg 3 
Mg 4 



- Mg +M 
+M 

- Mg 2 +M 
+M 

- Mg 3 +M 
+Mg 



Carbon 



1.20x10" 


32 





2160.0 


Same as 


F19 


4.40x10" 


10 





98600.0 


Same as 


F20 


4.40x10" 


10 





98600.0 


Same as 


F21 


2.76x10" 


29 





-2821.4 


Same as 


F22 


7.14x10 


-5 





17800 


Same as 


F23 


2.76x10" 


29 





-2821.4 


Same as 


F24 


1.66x10 


— 5 





19200 


Same as 


F25 


2.76x10" 


29 





-2821.4 


Same as 


F26 


8.30x10 


-7 





21600 


Same as 


F27 


8.30x10" 


10 








Same as 


F28 


8.30x10" 


10 








Same as 


F29 


- C 


4.36x10" 


18 


0.35 


161.3 


Andrcazza & Singh 1997 


1.00x10" 


17 








Clayton et al. 


1999 


1.00x10" 


10 








Clayton et al. 


1999 


1.00x10" 


13 








Clayton et al. 


1999 


1.00x10" 


10 








Clayton et al. 


1999 


1.00x10" 


10 








Clayton et al. 


1999 


1.00x10" 


13 








Clayton et al. 


1999 


1.00x10" 


10 








Clayton et al. 


1999 


1.00x10" 


■10 








Clayton et al. 


1999 


1.00x10" 


10 








Clayton et al. 


1999 


1.00x10" 


13 








Clayton et al. 


1999 


1.00x10" 


10 








Clayton et al. 


1999 


1.00x10" 


10 








Clayton et al. 


1999 


1.00x10" 


10 








Clayton et al. 


1999 


1.00x10" 


10 








Clayton et al. 


1999 


1.00x10" 


10 








Same as C2 


5.38x10" 


10 








Krusc & Roth 1997 


1.00x10" 


10 








Same as C3 


1.00x10" 


10 








Same as 


C18 


5.38x10" 


10 








Same as 


C17 


1.00x10" 


10 








Same as 


C18 


1.00x10" 


10 








Same as 


C18 


5.38x10" 


10 








Same as 


C17 


1.00x10" 


■10 








Same as 


C18 


1.00x10" 


10 








Same as 


C18 


5.38x10" 


10 








Same as 


C17 


1.00x10" 


10 








Curl & Haddon 1993 


1.00x10" 


10 








Same as 


C27 


1.00x10" 


10 








Same as 


C27 


1.00x10" 


10 








Same as 


C27 


1.00x10" 


■10 








Same as 


C27 


1.00x10" 


10 








Same as 


C27 


1.00x10" 


10 








Same as 


C27 


1.00x10" 


10 








Same as 


C27 



CI 

C2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

C10 

Cll 

C12 

C13 

C14 

C15 

C16 

C17 

C18 

C19 

C20 

C21 

C22 

C23 

C24 

C25 

C26 

C27 

C28 

C29 

C30 

C31 

C32 

C33 

C34 



C+ C 
C + C 2 
C + C 3 

c + c 4 



c 
c 
c 
c 
c 
c 



c 4 
c 5 
c 6 

C 6 

c 6 
c 7 



c + c 8 
c + c 8 
c 8 
c 8 
c 9 
■ c 2 

- c 2 

- c 3 

- c 4 
c 2 + c 4 
c 2 + c 5 

c 6 

C e 

c 7 
c 8 

C 8 

c 3 
c 4 

Cs 

c 6 
c 3 + c 7 
c 4 + c 4 
c 5 

C 6 



C - 
C - 
C - 

c 2 
c 2 
c 2 
c 2 



c 2 
c 2 
c 2 
c 2 
c 2 
c 3 
c 3 
c 3 
c 3 



c 4 
c 4 



C 2 + hi/ 
C 3 + hi/ 
C 4 + hi/ 



c 5 4 
c 2 4 
c 6 4 

C 7 4 
C 2 4 
C 3 -I 
C 8 4 
C 9 4 
C 2 4 
C 3 4 
C 4 4 
Cio 
C 4 4 
C 3 4 
C 5 4 
C 6 4 
C s 4 



hi/ 
C 3 
hi/ 
hi/ 
C 5 
C 4 
hi/ 
hi/ 
C 7 
C 6 
C 5 
- hu 
hi/ 
C 
hi/ 
hi/ 
C 



C 7 4- hi/ 
C 8 4- hu 



C 7 4 
C 9 4 
Cio 
C 9 4 
C 6 4 
C 7 4 
C 8 4 
C 9 4 
Cm 
C 8 4 
C 9 4 



C 
hi/ 

- hu 
C 
hi/ 
hi/ 
hi/ 
hi/ 

- hu 
hu 
hu 



Cio 4- hu 



Table 2 — Continued 



Chemical processes 



Aij 



Reference 



C35 
C36 
C37 
C38 
C39 
C40 
C41 
C42 
C43 
C44 
C45 
C46 
C47 
C48 
C49 
C50 
C51 
C52 
C53 
C54 
C55 
C56 
C57 
C58 
C59 
C60 



C5 + C5 
C+ C+ M 
C 2 + C 2 4 
O + C 2 
O + C 3 
O + C 4 

o + c 6 
o + c 6 
o + c 7 
o + c 8 
o + c 9 



M 



C 3 4 

C 4 4 

C 5 4 

C 6 4 

C 7 4 

C 8 4 

C 9 4 

C10 

He+ 

He+ 

He+ 

He+ 

He+ 

He+ 

He+ 



M 

M 

M 

M 

M 

M 

M 
r M 
4- C 3 
4- C 4 
+ C 5 
4- C 6 
4- C 7 
4- C 8 
4- C 9 



C10 + hu 
C 2 4- M 
C 4 4- M 
C 4 CO 
C 2 4- CO 
C 3 4- CO 
C 4 4- CO 
C 5 4- CO 
C 6 4- CO 
C 7 4- CO 
C 8 4- CO 
C 4 C 2 4M 
C 2 + C 2 +M 
C 2 4- C 3 4-M 
C 2 4- C 4 +M 
C 2 4- C 5 +M 
C 2 4- C 6 +M 
C 2 4- C 7 4-M 
C 2 4- C 8 +M 
C+ + C 2 4 He 
C+ 4- C 3 
C+ + C 4 
C+ 4- C 5 
C+ 4- C 6 
C+ 4- C 7 



He 
He 
He 
He 
He 



C+ 4- C 8 4 He 



l.OOxKT 10 Same as C27 

5.46X1CT 31 1.60 NIST 

2.48X1CT 21 -3.00 Kruse & Roth 1997 

5.99xl0~ 10 1420 NIST 

l.OOxlO -11 0.30 1130 Woon & Herbst 1996 

3.00X10- 10 0.00 1420 Tcrzieva & Herbst 1998 

l.OOxlO -11 -0.30 1130 Same as C39 

3.00X10- 10 0.00 1420 Same as C40 

l.OOxlO- 11 -0.30 1130 Same as C39 

3.00xl0~ 10 0.00 1420 Same as C40 

l.OOxlO -11 -0.30 1130 Same as C39 

6.64xl0~ 8 75406 Kruse & Roth 1997 

2.49xl0~ 8 69753 Reaction TF18 in CD09s 

6.64X10- 8 75406 Same as C46 

2.49xl0~ 8 69753 Same as C47 

6.64xl0~ 8 75406 Same as C46 

2.49xl0~ 8 69753 Same as C47 

6.64xl0~ 8 75406 Same as C46 

2.49xl0~ 9 69753 E 

1.60xl0~ 09 Reaction IM 15 in CD09 

1.60X10- 09 Same as C54 

1.60xl0-° 9 Same as C54 

1.60xl0~ 09 Same as C54 

1.60X10- 09 Same as C54 

1.60X10- 09 Same as C54 

1.60XHT 09 Same as C54 



a Thc rates arc given in the Arrhenius form kij(T) = Aij X (^g)" X cxp(— E a /T), where T is the gas 
temperature, Aij the Arrhenius coefficient in s _1 molecule -1 , cm 3 or cm 6 s _1 molecule -1 for a unimolecular, 
bimolccular or tcrmolccular processes respectively, v has no unit and reflects the temperature dependance of 
the reaction, and E a is the activation energy barrier in K _1 . 

b NIST is the NIST chemical kinetics database. UDFA06 is by Woodall ct al. (2007). Other references are 
listed in the bibliography. 'E' means 'estimated'. Reaction C53 has been estimated according to reaction C47 
decreased by a factor of 10 to account for the lower reactivity of the Cio ring compared to other chains. 

c The destruction rates by Compton electrons are calculated for the 170 Mq progenitor. For other progenitor 
masses, see Table 6 of CD09. 



Table 3. Gas temperature and number density profiles as a function of time for the 20 M & CCSN and the 170 M Q PISN used 
in this study. Similar profiles for the 20 M Q CCSN of Todini & Ferrrara (2001) and the 170 M PISN of Schneider et al. 

(2004) are also given for comparison. 







20 M 




170 M 


Todini & Fcrrara 


Schneider et al. 


Time (Days) 


T (K) 


iigas (cm~ 3 ) 


T (K) 


n gas (cm -3 ) 


T (K) 


rigos (cm 3 ) 


T(K) 


llgas (cm 3 ) 


100 


18000 


1.4 x 10 12 


21000 


3.9 x 10 11 


5400 


1.0 x 10 8 


5400 


1.0 X 10 12 


200 


5200 


1.8 x ll 11 


11800 


4.9 x 10 10 


3200 


1.0 x 10 8 


3200 


1.2 x 10 11 


300 


2500 


5.2 x 10 10 


6340 


1.4 x 10 10 


2400 


1.0 x 10 8 


2400 


3.7 x 10 10 


400 


1500 


2.2 x 10 10 


3330 


6.0 x 10 9 


1900 


1.0 x 10 8 


1900 


1.6 x 10 12 


500 


1000 


1.2 x 10 10 


1880 


3.1 x 10 9 


1600 


1.0 x 10 8 


1600 


8.0 x 10 9 


600 


740 


6.7 x 10 9 


1190 


1.8 x 10 9 


1400 


1.0 x 10 8 


1400 


4.6 x 10 9 


700 


560 


3.8 x 10 s 


857 


1.1 x 10 9 


1300 


1.0 x 10 8 


1300 


2.9 x 10 9 


800 


440 


2.8 x 10 9 


636 


7.5 x 10 8 


1100 


1.0 x 10 8 


1100 


1.9 x 10 9 


900 


360 


2.0 x 10 9 


479 


5.2 x 10 8 


1000 


1.0 x 10 8 


1000 


1.4 x 10 9 


1000 


300 


1.4 x 10 9 


470 


3.9 x 10 8 


960 


1.0 x 10 8 


960 


1.0 x 10 9 
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Table 4. Mass yields of clusters (in M ) ejected at day 1000 for the 170 and 20 M Q fully-mixed 

ejecta. a 



Clusters 


170 M 


20 M b 


(SiO) B 


25.39 


0.154 


(MgO) 4 








(FeO) 4 








AlO 


2.86 xl(T 3 


7.42 xl0~ 4 


(FcS) 4 








(MgS) 4 








(Fe) 4 


4.40xl(T 3 


4.32 xl0~ 4 


(Mg) 4 


3.18xl0~ 3 


1.56X10" 3 


(Si) 4 





1.59 xl0~ 4 


Cio 








(SiC) 2 








Total Cluster Mass 


25.4 


0.157 


Efficiency 


30.9 % 


2.7 % 



a The efficiency is denned as the ratio of the cluster 
mass to the He core mass given in CD09. 

b A mass cut of 2.4 Mq is assumed as in Nozawa 
ct al. (2003). 
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Table 5. Dust mass yields (in M ) ejected at day 1000 for the 170 M Q and the 20 M Q 
fully-mixed ejecta. The dust mass yields from Nozawa et al. (2003) for the same mass progenitors 
are also listed. Values quoted for the present study are upper limits to dust yields. Two cases of 
metal depletion are considered: a case where Mg/Fe atoms are totally depleted by dust formation 
and a case where no metal depletion occurs (see text). The % columns represent the mass 
fraction of the total dust mass for each ejecta. 



This Study NK03 



Dust 


Mg/Fe depleted 


% of Dust 


No depletion 


% of Dust 


Dust Mass 


% of Dust 


170 M 














Silica (Si0 2 ) 


12.9 


38.7 % 


17.3 


68.1 % 


20.0 


56.1 % 


Fayalitc (Fc 2 Si0 4 ) 


6.6 


19.8 % 












Forsterite (Mg 2 Si04 ) 


5.7 


17.1 % 







11.3 


31.7 % 


Pure silicon (Si) 


8.1 


24.3 % 


8.1 


31.8 % 







Alumina (AI2O3) 


0.003 


0.009 % 


0.003 


0.012 % 


0.05 


0.14 % 


Magnetite (FesCU) 












4.3 


12.0 % 


Pure iron (Fe) 


0.004 


0.012 % 


0.004 


0.016 % 







Carbon (AC) 

















Silicon carbide (SiC) 

















Total Dust 


33.3 


100 % 


25.4 


100 % 


35.6 


100 % 


20 M 














Silica (Si0 2 ) 







0.105 


67.5 % 


0.42 


57.5% 


Fayalite (Fe 2 Si0 4 ) 


0.125 


37.4 % 












Forsterite (Mg 2 Si0 4 ) 


0.160 


47.7 % 







0.23 


31.5 % 


Pure silicon (Si) 


0.049 


14.7 % 


0.049 


31.5 % 







Alumina (AI2O3) 


8.8 xl0~ 4 


0.26 % 


8.8 xl0~ 4 


0.56 % 


9.0 xl0~ 4 


0.12 % 


Magnetite (Fe304) 












0.08 


10.9 % 


Pure iron (Fe) 







4.32 xl0~ 4 


0.28 % 







Pure Magnesium (Mg) 







0.0014 


0.88 % 







Carbon (AC) 

















Silicon carbide (SiC) 

















Total Dust 


0.33 


100 % 


0.16 


100 % 


0.73 


100 % 



Table 6. Molecular and dust yields (in M Q ) ejected at 1000 days for our comparative study of Schneider at al. (2004, SFS04) 
models. Values quoted for the present study are upper limits to dust yields. Cluster yields are listed under 'Ejecta". For each 
progenitor, two cases of metal depletion have been considered (see text) : 1) Mg and Fe are totally depleted by dust formation, 

and 2) no metal depletion is considered. 







170 M Q 








260 M 






Molecules / Clusters /Dust 


Ejecta 


Mg/Fe depleted 


No depiction 


SFS04 


Ejecta 


Mg/Fe depleted 


No depletion 


SFS04 


CO 


10.5 


10.5 


10.5 


5 


8.5 


8.5 


8.5 


0.01 


SO 


9.7 


9.7 


9.7 




17.8 


17.8 


17.8 






30.1 


13.3 


30.1 




9.6 


6.3 


9.6 




lotai (ivioiccuicsj 


OU.O 


oo.O 


OU.O 


a 



qc n 

oo.y 


QO ft 


oo.y 


U.U1 




z / .0 








QCl P. 

oy .0 












































a i r\ 


Z.oZ X 1U 








Z.Zz X 1U 


























(MgS) 4 


_ 








_ 








(Fe) 4 


7.28xl0~ 4 








34.3 








(Mg) 4 


1.78X10- 1 








2.5 








(Si) 4 


















Cio 


















Total (Clusters) 


27.8 








74.8 








Silica (Si0 2 ) 




12.1 


18.8 


26 




14.4 


25.9 


46.0 


Fayalite (Fe 2 Si0 4 ) 




4.7 








12.6 






Forsterite (Mg 2 Si0 4 ) 




12.3 




11.0 




6.0 




11.0 


Pure silicon (Si) 




8.8 


8.8 






12.1 


12.1 




Alumina (AI2O3) 




0.03 


0.03 


0.03 




0.02 


0.02 


0.03 


Magnetite (Fe30 4 ) 








0.06 








8.7 


Pure iron (Fe) 







7.28xl0~ 4 






34.3 


34.3 




Pure iron (Mg) 







1.8 x 10- 1 






2.5 


2.5 




Carbon (AC) 








1.6 








3.3 


Total (Dust) 




37.9 


27.8 


38.7 




81.9 


74.8 


69.1 
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Table 7. Upper limits to dust mass yields (in M & ) ejected at day 1000 for the 170 M unmixed ejecta. Two cases of metal 
depletion are considered: 1) a case where Mg is totally depleted by dust formation (D), and 2) a case where no metal depletion 

occurs (ND) - see text for more detail. a 





Zone 1 


Zone 2 


Zone 3 


Zone 4 




Zone 5 




Zone 1-5 


Zone mass 


(20 M e ) 


(15 Mo) 


(15 


M ) 


(23 M ) 




(4.3 Mas) 




(82.3 M ) 


Major elements' 3 


Si/S/Fe 


O/Si/S/Mg 


O/Mg/Si 


O/C/Mg 




o/c 










ND 


D 


ND 


D 


ND 


D 


ND 


D 


ND 




D 


ND 


Silica 




2.577 


2.976 





0.661 





1.8xl0~ 3 





3.2x10" 


-6 


2.577 


3.638 


Forsterite 




0.927 





1.542 





4.3xl0~ 3 





7.5x10 


-6 




2.474 





Alumina 




9.2xl0~ 4 


9.2X10- 4 


0.028 


0.028 


8.0xl0- 4 


8.0xl0~ 4 










0.0296 


0.0297 


Pcriclasc 




2.8xl0~ 9 


2.8xl0~ 9 


2.5xl0~ 6 


2.5xl0~ 6 
















2.5xl0~ 6 


2.5xl0~ 6 


Pure iron 


6.7xl0~ 5 




















6.7xl0~ 5 


6.7xl0~ 5 


Pure silicon 


0.265 


1.388 


1.388 


0.308 


0.308 


8.6xl0~ 4 


8.6xl0~ 4 


1.5x10 


~ 6 1.5x10" 


-6 


1.963 


1.963 


Pure magnesium 




3.1xl0~ 4 


3.1xl0~ 4 


7.4xl0~ 3 


7.4xl0~ 3 


2.5xl0~ 4 


2.5xl0~ 4 


3.3x10 


" 8 3.3x10" 


-8 


2.5xl0~ 4 


0.00802 


Iron sulphide 


0.011 




















0.011 


0.011 


Carbon 












1.5xl0~ 14 


1.5xl0~ 14 


1.1x10" 


" 13 1.1x10" 


13 


1.2xl0~ 13 


1.2X10 -13 


Silicon carbide 


























Total Dust 


0.277 


4.894 


4.365 


1.885 


1.004 


0.0062 


0.0038 


9.0x10 


~ 6 4.7x10" 


-6 


7.06 


5.65 


Efficiency 


1.38 % 


32.63 % 


25.1 % 


12.57 % 


6.69 % 


0.027 % 


0.016 % 


% 


% 




8.58 % 


6.86 % 



a The efficiency is defined as the ratio of the molecular mass to the zone mass. 
b In order of decreasing mass yield (see Table 3 of CD09). 



Table 8. Upper limits to dust mass yields (in M & ) ejected at day 1000 for the 20 M & unmixed ejecta. Two cases of metal 
depletion are considered: 1) a case where Mg is totally depleted by dust formation (D), and 2) a case where no metal depletion 

occurs (ND) - see text for more detail. a 





Zone 1 


Zone 2 


Zone 3 




Zone 4 


Zone 1-4 


Zone mass 


(0.6 Ma,) b 


(0.6 M Q ) 


(1.35 M©) 




(0.9 Mr,) 


(3.45 M Q ) 


Major elemcnts c 


Si/S/Fe 


O/Mg/Si 


O/C/Mg 




Hc/C/O 








ND 


D 


ND 


D 


ND 


D 


ND 


D 


ND 


Silica 


1.4xl(r 3 





0.038 





2.8 xl0~ 5 











0.039 


Forsterite 




0.089 





6.5 xl0~ 5 











0.089 





Alumina 




7.9xl0~ 5 


7.9xl0~ 5 


7.4xl0- 6 


7.4xl0- 6 








8.6xl0~ 5 


8.6xl0~ 5 


Periclase 




4.9xl0~ 4 


4.9xl0~ 4 














4.9xl0~ 4 


4.9xl0~ 4 


Pure iron 


4.6X1CT 5 














4.6xl0~ 5 


4.6xl0~ 5 


Pure silicon 


0.012 


0.018 


0.018 


1.3X10" 6 


1.3xl0" s 








0.030 


0.030 


Pure magnesium 




3.9xl0~ 4 


3.9xl0~ 4 


2.1xl0~ 6 


2.1xl0~ 6 








3.9xl0~ 4 


3.9xl0~ 4 


Iron sulphide 


0.033 














0.033 


0.033 


Carbon 








1.8X10 -15 


1.8X10- 15 


1.3x10 


" 13 1.3xl0~ 13 


1.3xl0- 13 


1.3X10- 13 


Silicon carbide 




















Total Dust 


0.046 


0.108 


0.057 


8.8xl0 -5 


5.1xl0" 5 








0.154 


0.103 


Efficiency 


7.67 % 


18.0 % 


9.5 % 


0.0065 % 


0.0038 % 


% 


0% 


4.47 % 


2.98 % 



a The efficiency is denned as the ratio of the molecular mass to the zone mass. 
b A mass cut of 2.4 Mq is assumed for Zone 1 as in NK03 
c In order of decreasing mass yield (see Table 3 of CD09). 
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Table 9. Depletion efficiencies of metals in molecular clusters and molecules (underlined) for the 

170 and 20 M unmixed ejecta. a 



170 M 20 M 



Element 


Clusters 


Zone 1 


Zone 2 


Zone 3 


Zone 4 


Zone 5 


Zone 1 


Zone 2 


Zone 3 


Zone 4 


Si 


(Si0 2 )s 




50.0 % 


50.0 % 


46.7 % 


37.2 % 


0.12 % 


49.5 % 


48.2 % 






(Si) 4 ,5 


1.8 % 


50.0 % 


50.0 % 


46.7 % 


37.2 % 


3.1 % 


49.5 % 


48.2 % 






SiS 


26.5 % 










28.3 % 








O 


(Si0 2 ) 5 




9.6 % 


2.7 % 








3.9 % 








o 2 




38.0 % 


46.4 % 


12.3 % 






60.9 % 


27.6 % 






CO 








29.1 % 


43.0 % 






32.7 % 






so 




L2 % 


03 % 


Ol % 




07% 








s 


(FcS) 4 


0.09 % 










8.7 % 










SiS 


98.9 % 










91.4 % 










SO 




100 % 


99.8 % 








99.2 % 






Al 


AlO 




97.0 % 


98.3% 


92.1% 






99.1 % 


96.8 % 




Fe 


(FeS) 4 


1.8 % 










44.0 % 










(Fe) 4 


0.02 % 










0.1 % 








Mg 


(Mg) 4 




0.1 % 


2.3 % 


0.1 % 






0.91 % 


0.04 % 






(MgO) 4 
















0.73 % 




C 


Cio 










- 21.6 % b 








- 95.3 % b 




CO 






100.0 % 


100.0 % 


77.3 % 




99.8 % 


99.7 % 


- 0.14 % b 



a The depletion efficiency is defined by Eq. (41) in the text 
b The efficiencies are for a He-rich and He-free Zone. 
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Fig. 1. — The structures of some small clusters are illustrated for silicone monoxide (adapted from 
Lu et al. 2006), magnesium monoxide (adapted from Kohler et al. 1997), and pure magnesium 
(adapted from Jellinek & Acioli 2002). Structure similar to (MgO) n applies to iron oxide (FeO) n 
and magnesium sulphide (MgS) n and iron sulphides (FeS) n . Pure silicon and iron clusters have 
structures similar to pure magnesium (Mg) n clusters. 
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Fig. 2. — The evolution of the cluster abundances normalized to total gas number density for the 
170 Mq fully-mixed ejecta: a) Silicon oxide clusters and aluminum monoxide; b) Iron clusters, most 
important ions, and electrons. The dominant precursors are silicon monoxide, pure magnesium and 
iron clusters, and molecular aluminum monoxide. Carbon chains and rings do not form in fully- 
mixed ejecta. The abundance of secondary electrons is commanded by the ionization of the metals 
Mg and Fe. The He + abundance is high from 100 to 400 days impeding the efficient synthesis of 
molecules during this timespan (see CD09). The O2 abundance is also shown. 




Fig. 3. — The evolution of the cluster abundances normalized to total gas number density for the 
20 Mq fully-mixed ejecta: a) Silicon oxide clusters and aluminum monoxide; b) Iron clusters, most 
important ions, and electrons. Carbon chains/rings do not form in significant amounts. The He + 
abundance is high from 100 to 700 days impeding the efficient synthesis of molecules during this 
timespan (see CD09). The abundance of secondary electrons is commanded by the ionization of 
He and metals like Mg and Fe. The O2 abundance is also shown. 
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Fig. 4. — Carbon clusters abundances normalized to total gas number density as a function of 
post-explosion time for the 20 M & fully-mixed ejecta model of Todini & Ferrara (2001, TF01). 
Carbon chains abundances (with respect to total gas density) are extremely small due to the large 
amounts of He + in the ejecta. These results totally contradict the findings by TF01 who predict 
the formation of amorphous carbon grains in the ejecta of primordial SNe with 20 M & progenitors 
- see §4.2.2. 
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Fig. 5. — The evolution of the cluster abundances normalized to total gas number density for the 
170 Mq unmixed ejecta - Zone 1: a) Silicon clusters - molecular SiS and SiO are also shown; b) 
Iron-bearing clusters, most important ions, and electrons. The dominant precursors are pure silicon 
and iron sulphide clusters. The abundance of secondary electrons is commanded by the ionization 
of atomic Fe. 




Fig. 6. — The evolution of the cluster abundances normalized to total gas number density for 
the 170 Mq unmixed ejecta - Zone 2: a) Silicon monoxide clusters; b) Magnesium clusters and 
aluminum monoxide, most important ions, and electrons. The dominant precursors are silicon 
monoxide clusters and aluminum monoxide. The abundance of secondary electrons is commanded 
by the ionization of atomic Mg. 
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Fig. 7. — Dust chemical composition and mass yields (in units of Mq) at t = 1000 days versus 
depletion case for the 170 Mq unmixed ejecta. Results from Nozawa et al. (2003) are shown for 
comparison. The total upper limit for dust mass and the total dust mass from NK03 are also 
indicated. Zone 1 synthesizes pure Si, FeS, and Fe clusters. Silica and/or silicates (forsterite), pure 
silicon clusters are formed in Zone 2 and Zone 3. No dust is formed in Zone 4. Very small amounts 
of silica and/or silicates (forsterite), and pure silicon clusters form in Zone 5, when AC and SiC 
grains do not. 
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Fig. 8. — The evolution of the cluster abundances normalized to total gas number density for the 
20 Mq unmixed ejecta - Zone 1: a) Silicon-bearing clusters - Si/S-bearing species; b) Iron-bearing 
clusters, most dominant atoms and ions, and electrons. The dominant precursors are pure silicon, 
iron sulphide, silicon monoxide and pure iron clusters. The abundance of secondary electrons is 
commanded by the ionization of atomic Fe. 




Fig. 9. — The evolution of the cluster abundances normalized to total gas number density for 
the 20 Mq unmixed ejecta - Zone 2: a) Silicon monoxide clusters, aluminum-bearing species; b) 
Magnesium-bearing clusters, atoms, ions and electrons. The dominant precursors are silicon oxide 
and pure magnesium clusters, and aluminum monoxide. The abundance of secondary electrons is 
commanded by the ionization of atomic Mg. 
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Fig. 10. — Dust chemical composition and mass yields (in units of Mq) at t = 1000 days versus 
depletion case for the 20 Mq unmixed ejecta. Results from Nozawa et al. (2003) are shown for 
comparison. The total upper limit for dust mass and the total dust mass from NK03 are also 
indicated. Zone 1 synthesizes iron sulphide, pure Si clusters, and Fe clusters. Silica or forsterite, 
pure silicon clusters and alumina are formed in Zone 2 and Zone 3. Extremely small amounts of 
carbon chains are formed in Zone 4. 




Fig. 11. — Abundances normalized to total gas number density of (SiO)s clusters (solid lines) and 
gas-phase SiO (dotted lines) as a function of time in Zone 2 of the 170 M@ unmixed ejecta. Curves 
are labelled by the decrease factor applied to reaction [7] rates according to similar decrease in total 
gas pressures: 1 corresponds to the rate at 1 atm, while results presented in §4.2 correspond to a 
decrease of 10~ 2 . For decreasing nucleation rates, the formation of (SiO)s clusters is delayed to 
late times though the final mass yield for the clusters is similar for all cases. As expected, the SiO 
abundance is larger and its decreases with time is less pronounced for the smaller reaction [7] rate. 
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Fig. 12. — Abundances for Cio rings normalized to total gas number density for various He initial 
composition in Zone 4 of the 20 M & unmixed ejecta. Curves are labelled by the percentage of the 
initial He mass present. The 100 % labelled curve corresponds to the results for Zone 4 discussed 
in §4.3.2. The larger the He content, the lower the abundance of the end-product ring Cio due to 
the destruction of C-chains by He + , and the larger the delay in formation time, which depends on 
the He + abundance in the ejecta. 



